•  • 

NUWC’NPT  Technical  Raport  10.539 
23  Novambar  1993 


AD-A273  999 

■ililiiii 

Detection  Performance  of  a  Modified 
Generaiized  Likeiihood  Ratio 
Processor  for  Random  Signals  of 
Unknown  Location 


Albert  H.  Nuttall 
Surface  A8W  Directorate 


in' 

i  I V.* 

DEC 2  2  1^93  "  l\ 
'  k  i 


Naval  U  ndersea  Warfare  Center  Detachment 

New  London,  Connecticut 


Approved  for  public  release;  distribution  Is  unlimited. 


93-3084 1 

Ill/ll  fllll  llfll  Hill  aain  ■ 


93  12  2l  I8l 


PRBFACB 


This  research  was  conducted  under  NUHC  Job  Order  Nuad>er 
A10020,  R&D  Project  Munber  RROONOO,  Perforaance  Evaluation  of 
EOnlinear  Signal  Processors  with  Misaatch,  Principal  Investigator 
Dr.  Albert  H.  Nuttall  (Code  302).  This  technical  report  was 
prepared  with  funds  provided  by  the  NUWC  In-House  Independent 
Research  Progran,  sponsored  by  the  Office  of  Naval  Research. 

Also,  the  research  presented  in  this  report  was  sponsored  by 
the  Technology  Directorate  of  the  Office  of  Naval  Research, 

T.  G.  (k>ldsberry  (ONR  451). 


The  author  of  this  report  is  located  at  the  Naval  Undersea 
Warfare  Center  Detachment,  New  London,  CT  06320.  The  technical 
reviewer  for  this  report  was  Erilc  B.  Siborg,  (Code  2111). 


REVIEWED  AND  APPROVED:  23  Roveid)er  1993 


Donald  W.  Counsellor 
Director,  Surface  Antisubmarine  Warfare 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No.  0704-0188 

*«alK  wBonim  auMtn  *Of  Uwt  comttiew  e(  iMormiuon  <\  mimaifd  te  <  xouf  ott  rttpomc.  includHtg  tlw  imc  ter  raviewHiq  untructiom.  tMrcfunq  uniuif  deu  Murcm, 

jne«fiii9  iM  nnwiiiWMin  itie  eew  wMded.  and  teweneeq  «««■  r«»m»im  tiw  tawtctieii  el  inlornuiien  tettdcommtftareeardrnf  iluiburdtn  ttmuner  wiy  oilitr  moki  e<  itw 
cemaiew  et  HilenMnen.  mcluduiQ  tuMtiuem  tor  reduonq  ten  Burd«i.  lo  WMiunqtOA  Mdaeoudrifn  tervK«.  OrrcctetaM  ter  mtermauen  Oeeranent  and  Ueeru.  )21S  iattenon 
OawiHifiMiav.Vwie  in*.  &rlin9ian.V»  1220I-4N2.  andieiKdOttKtetManaddmantandtuddat.  daearwarktadu(iienrreiect(0>04-aiW).Wailwi9ien.  DC  20M}. 

1.  AGiNCY  USf  ONLY  a««v*  blank)  2.  RiPORT  OATi  3.  REPORT  TYPE  AND  OATES  COVERED 

23  November  1993  Progress 

4.  TITU  AND  SUtTITLE 

Detection  Performance  of  a  Modified 

Generalized  Likelihood  Ratio  Processor 
for  Random  Signals  of  Unknovm  location 

S.  FUNDING  NUMBERS 

PE  0601152N 

C.  AUTHOR(S) 

Albert  H.  Nuttall 

7.  PERFORMING  ORGANUATION  NAME(S)  AND  AOORESS(ES) 

Naval  Undersea  Warfare  Center  Detachment 

New  London,  Connecticut  06320 

B.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

NUWC-NPT  TR  10,539 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AODRESS<ES) 

Office  of  Naval  Research 

800  North  Quincy  Street,  BCT  1 

Arlington,  VA  22217-5000 

10.  SPONSORING /MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

12a.  DISTRIBUTION /AVAILA8IUTY  STATEMENT 

Approved  for  public  release; 
distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (MaMimum  200  words) 

A  search  region  of  N  bins  in  time,  frequency  space  contains 
either  noise  only,  or  noise  plus  a  signal  which  occupies  a  subset  of 

M  bins  of  unknown  location.  The  unknown  signal  strengths  are 
estimated  by  means  of  a  maximum  likelihood  procedure  modified  by  a 
breakpoint  x  ,  in  order  to  suppress  the  weaker  estimates.  The 
stronger  signal  estimates  are  then  used  to  implement  a  modified  form 
of  generalized  likelihood  ratio  processor. 

The  highly  nonlinear  character  of  this  generalized  likelihood 
ratio  processor  debilitates  against  an  exact  analysis  of  its 
perfozmance.  To  circumvent  this  difficulty,  all  the  nonlinear 
components  are  replaced  by  piecewise-linear  devices  which  best  match 
the  given  characteristics  according  to  minimum  mean  square  error. 

The  false  alarm  and  detection  probabilities  of  this  piecewise-linear 
processor  are  then  derived  exactly.  These  analytic  results  enable 
ready  determination  of  threshold  settings  and  receiver  operating 

'(le1:ect?on  performance  modified  processor 

generalized  likelihood  ratio 

random  signals  unknown  location 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

UNCLASSIFIED 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

SAR 

MSN  7540-01 -280-5500  Stanoard  form  298  (9ev  2-89) 


^'•%CriOM  Ov  ^NSi  >ld  ^39*  S 


:99-i03 


UHCXASSZFIBD 


OF  THIS  PAGE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION 
OF  THIS  PAGE 


TR  10539 


TABLE  OF  CONTENTS 


Page 


LIST  OF  FIGURES  ill 

LIST  OF  TABLES  V 

LIST  OF  SYMBOLS  vi 

INTRODUCTION  1 

PROBLEM  DEFINITION  3 

DERIVATION  OF  FALSE  ALARM  PROBABILITY  5 

Approach  5 

Piecewise-Linear  J^proximation  5 

Exceedance  Probability  for  Noise  Alone  7 

DERIVATION  OF  DETECTION  PROBABILITY  11 

Statistics  of  Bin  Outputs  with  Signal  11 

Statistics  of  Summer  Output  C  12 

FOURTH-ORDER  MOMENT  APPROACH  17 

ENERGY  DETECTOR  PERFORMANCE  21 

GRAPHICAL  RESULTS  23 

Piecewise-Linear  Fit  h(x)  23 

False  Alarm  Probability  26 

Detection  Probability  33 

Variation  of  Performance  with  N  73 

Energy  Detector  Characteristics  81 

EXTENSION  TO  UNEQUAL  SIGNAL  PONERS  91 

Matching  of  Mean  Outputs  92 

Matching  of  First  Two  Moments  94 


SUMMARY 


101 


TR  10539 


APPENDICES 

A.  MODIFIED  GENERALIZED  LIKELIHOOD  RATIO  PROCESSOR  103 

B.  PIECEWISE-LINEAR  FIT  TO  NONLINEARITY  g(x)  107 

C.  PARTIAL  FRACTION  EXPANSION  OF  PRODUCT  117 

D.  TWO  ALTERNATIVE  EXACT  ANALYSIS  METHODS  119 

High-Order  Moment  Approach  119 

Characteristic  Function  Approach  120 

Asymptotic  Behavior  of  Characteristic  Function  121 

Evaluation  of  Derivatives  of  Py(u)  124 

E.  TABLES  OF  FIRST  FOUR  CUMULANTS  OF  y.  127 

'n 

F.  INDEFINITE  INTEGRAL  OF  EXPONENTIAL  INTEGRAL  131 

G.  DIRECT  EVALUATION  OF  EXACT  EXCEEDANCE  DISTRIBUTION  133 

REFERENCES  135 


ii 


TR  10539 


LIST  OF  FIGURES 


Figure 

1 .  FIT-2  for 

S 

-  0  - 

Page 

24 

2. 

3. 

FIT-2  for 

?£  versus 

*o  "  5 

V  for  N 

a 

1024 

,  FIT-2 

25 

28 

4. 

?£  versus 

V 

for  M 

- 

512, 

FIT-2 

29 

5. 

versus 

V 

for  N 

- 

256, 

FIT-2 

30 

6. 

?£  versus 

V 

for  N 

- 

128, 

FTT-2 

31 

7. 

?£  versus 

V 

for  M 

- 

64, 

FIT-2 

32 

8. 

ROC  for  N 

• 

1024, 

N 

■  8, 

-  1,  FIT-4 
o 

37 

9. 

ROC  for  N 

m 

1024, 

N 

■  8# 

Xq  -  3,  FIT-4 

38 

10. 

ROC  for  N 

m 

1024, 

M 

■  8, 

•  5,  FIT-2 

39 

11. 

ROC  for  N 

m 

1024, 

M 

-  8, 

X^  -  7,  FIT-2 

40 

12. 

ROC  for  N 

m 

1024, 

M 

•  8, 

X^  •  9,  FIT-2 

41 

13. 

ROC  for  M 

m 

1024, 

N 

■  8, 

Xj,  -  11,  FIT-2 

42 

14. 

ROC  for  N 

m 

1024, 

M 

■  16 

,  Xq  -  1,  FIT-4 

43 

15. 

ROC  for  N 

m 

1024, 

N 

■  16 

,  x„  -  3,  FIT-4 
o 

44 

16. 

ROC  for  N 

m 

1024, 

M 

■  16 

,  x^  -  5,  FIT-2 
o 

45 

17. 

ROC  for  N 

m 

1024, 

N 

-  16 

,  Xq  -  7,  FIT-2 

46 

18. 

ROC  for  N 

m 

1024, 

N 

■  16 

,  Xq  -  9,  FIT-2 

47 

19. 

ROC  for  N 

m 

1024, 

M 

-  16 

,  Xjj  -  11,  FIT-2 

48 

20. 

ROC  for  N 

m 

1024, 

M 

-  32 

,  Xq  -  1,  FIT-4 

49 

21. 

ROC  for  N 

m 

1024, 

M 

-  32 

,  x„  -  3,  FIT-4 

50 

22. 

ROC  for  N 

m 

1024, 

M 

«  32 

,  Xq  -  5,  FIT-2 

51 

23. 

ROC  for  N 

m 

1024, 

M 

«  32 

,  Xq  -  7,  FIT-2 

52 

24. 

ROC  for  N 

m 

1024, 

M 

-  32 

,  Xp  -  9,  FIT-2 

53 

•  •  • 
1X1 


TR  10539 


Figure 

Page 

25. 

ROC 

for 

N 

«l 

1024,  N 

-  32, 

Xq  -  11,  FIT-2 

54 

26. 

ROC 

for 

N 

m 

1024,  M 

-  64, 

x^  -  1,  FIT-4 

55 

27. 

ROC 

for 

N 

m 

1024,  M 

-  64, 

Xj,  -  3,  FIT-4 

56 

28. 

ROC 

for 

N 

m 

1024,  M 

-  64, 

Xq  -  5,  FIT-4 

57 

29. 

ROC 

for 

N 

m 

1024,  M 

-  64, 

Xq  -  7,  FIT-2 

58 

30. 

ROC 

for 

N 

m 

1024,  M 

-  64, 

Xq  -  9,  FIT-2 

59 

31. 

ROC 

for 

N 

m 

1024,  M 

-  64, 

Xq  -  11,  FIT-2 

60 

32. 

ROC 

for 

N 

m 

1024,  M 

-  128, 

Xjj  -  1,  FIT-4 

61 

33. 

ROC 

for 

N 

m 

1024,  M 

-  128, 

Xjj  -  3,  FIT-4 

62 

34. 

ROC 

for 

N 

m 

1024,  M 

-  128, 

Xjj  -  5,  FIT-4 

63 

35. 

ROC 

for 

N 

m 

1024,  M 

-  128, 

Xq  -  7,  FIT-2 

64 

36. 

ROC 

for 

N 

m 

1024,  M 

«  128, 

X^  -  9,  FIT-2 

65 

37. 

ROC 

for 

N 

m 

1024,  M 

»  128, 

Xq  •  11,  FIT-2 

66 

38. 

ROC 

for 

N 

m 

1024,  N 

-  256, 

Xq  -  1,  FIT-4 

67 

39. 

ROC 

for 

N 

m 

1024,  M 

-  256, 

Xq  -  3,  FIT-4 

68 

40. 

ROC 

for 

N 

m 

1024,  N 

-  256, 

Xq  -  5,  FIT-4 

69 

41. 

ROC 

for 

N 

m 

1024,  M 

«  256, 

Xq  -  7,  FIT-2 

70 

42. 

ROC 

for 

N 

m 

1024,  M 

-  256, 

Xq  -  9,  FIT-2 

71 

43. 

ROC 

for 

N 

m 

1024,  N 

-  256, 

-  11,  FIT-2 

72 

44. 

ROC 

for 

M 

m 

512,  N  - 

*0 

-  11,  FIT-2 

75 

45. 

ROC 

for 

N 

m 

256,  N  - 

*o 

■  11,  FIT-2 

76 

46. 

ROC 

for 

N 

m 

512,  N  - 

32,  X 

^  -  7,  FIT-2 
o 

77 

47. 

ROC 

for 

N 

- 

256,  M  - 

32,  z 

^  •  7,  riT-2 
o 

78 

48. 

ROC 

for 

N 

- 

512,  M  - 

256, 

Xq  -  1,  FIT-4 

79 

49. 

ROC 

for 

N 

- 

256,  N  - 

256, 

Xq  -  1,  FIT-4 

80 

IV 


TR  10539 


Figure  Page 

50.  Energy  Detector  False  Alarsi  Probability  82 

51.  Energy  Detector  ROC  for  N  •  1024,  N  •  8  83 

52.  Energy  Detector  ROC  for  M  •  1024,  M  ■  16  84 

53.  Energy  Detector  ROC  for  N  «  1024,  N  «  32  85 

54.  Energy  Detector  ROC  for  N  •  1024,  M  ■  64  86 

55.  Energy  Detector  ROC  for  N  ■  1024,  M  >  128  87 

56.  Energy  Detector  ROC  for  M  ■  1024,  M  *  256  88 

57.  Energy  Detector  ROC  for  M  •  1024,  M  ■  512  89 

58.  Energy  Detector  ROC  for  N  >  1024,  N  ■  1024  90 

59.  Hean  versus  Signal  Power  S  97 

60.  Difference  of  Means  versus  Signal  Power  S  98 

61.  Difference  of  Variances  versus  Signal  Power  S  99 

62.  Ratio  R(S,Xq)  versus  Signal  Power  S  100 

G-1  False  Alarm  Probability  for  N  >  1024,  ■  1  134 


LIST  OF  TABLES 

Table  1.  S(dB)  Required  for  N  ■  1024,  P^  ■  10”^,  P^j  ■  .5  36 
Table  2.  S(dB)  Required  for  N  •  1024,  P^  ■  10”®,  Pj|  ■  36 
Table  3.  S(dB)  Required  for  P^  -  10”^,  P^  -  .5  74 
Table  4.  S(dB)  Required  for  P^  ■  10”®,  P^j  ■  *8  74 


V 


TR  10539 


M6LR 

N 


n 

bold 


“0 

^0 


Eixnl 


overbar 
M 
H, 


‘1 


o 

g(x) 


'n 


z 

V 

h(x) 

c,b 

C 

Px 


LIST  OF  SYMBOLS 

modified  generalized  likelihood  ratio 
total  number  of  search  bins,  (1) 
n-th  bin  measurement  or  observation,  (1) 
random  variable 
noise-only  hypothesis,  (1) 

probability  density  function  of  under  Hq,  (1) 
ensemble  average  of  z^,  (1) 
ensemble  average,  (1),  (44) 

number  of  bins  occupied  by  signal  (when  present),  (2) 

signal -plus-noise  hypothesis,  (2) 

probability  density  function  of  z^^  under  Hj,  (2) 

average  signal  power  in  n-th  bin,  (2) 

signal  power  parameter  (1  +  S^^)"^,  (2) 

minimum  acceptable  signal  power  level  estimate,  (3) 

breakpoint  1  S^,  (3) 

nonlinear  function,  (3) 

n-th  distorted  random  variable  g(*,j)#  (4) 
decision  variable,  processor  output,  (4) 
threshold,  (4) 

approzimating  nonlinearity,  (5) 

constants  in  device  h(z),  (5) 

output  of  device,  h(Zjj),  (7) 

output  of  summer,  approximation  to  z,  (8) 

false  alarm  probability 

detection  probability 

vi 


TR  10539 


o  o 

Pr(B) 

B 

Pn 

U(x) 

fw 


P 


n 


n 


Cn 
Pc 

=C 

^C 

Qi 
s 

a 

A 

K 

f.pfi 


‘a'‘b 


’’k 

^kn 

'’mn 


parameters  of  device  h(x)  for  noise  only,  (10) 
exceedance  distribution  function  of  (10) 
probability  of  event  B,  (10) 
exp(-x^),  (11) 

probability  density  function  of  (12) 
unit  step:  1  for  x  >  0,  0  otherwise,  (12) 
argument  of  characteristic  function 
characteristic  function  of  device  output  >|^,  (13) 
characteristic  function  of  summer  output  Cr  (14) 
normalised  probability  density  function,  (15) 
normalized  characteristic  function,  (16) 
normalized  exceedance  distribution  function,  (17) 
normalized  cumulative  distribution  function,  (18) 
probability  density  function  of  C,,  (19) 
exceedance  distribution  function  of  C#  (20) 

cumulative  distribution  function  of  (21) 

exceedance  distribution  function  of  x^  under  H^,  (22) 
equal  signal  power  level,  (22) 
signal  power  parameter,  1/(1  S),  (22) 

exp(-axQ),  (24) 

N  -  M,  (27) 

component  characteristic  functions,  (27) 
expansion  coefficients,  (29) 
expansion  coefficients,  (31) 
expansion  coefficients,  (33),  (37) 
expansion  coefficients,  (33),  (37) 


TR  10539 


Xy(lt) 

X,(k) 
Py(J«) 
f4(5) 
b'  ,c' 
w'  ,M' 
X,. 


w 

ROC 

ED 

PCS/Xq) 

E^(z) 

% 

R(S,X„) 


z 

FIT-2 

FIT-1 

*z<S' 


k-th  cufflulant  of  y^,  (43) 
k-th  cumulant  of  z,  (43) 
k-th  moment  of  (44) 

candidate  characteristic  function,  FIT-4,  (46) 

numerator  constants  in  f^((),  (46) 

denominator  constants  in  f^(C),  (46) 

auxiliazy  constants,  (47) 

energy  detector  output,  (52) 

characteristic  function  of  w,  (53) 

receiver  operating  characteristic,  figure  8 

maximum  false  alarm  probability,  (59) 

energy  detector,  table  1 

mean  output  of  g(x),  (63) 

exponential  integral,  (63) 

mean  of  summer  output  x,  (65) 

effective  signal  power,  (66),  (67) 

effective  number  of  signal  bins,  (68) 

peak  signal  power,  max{^),  (68) 

variance  of  output  of  g(x),  (69) 

ratio  of  functions,  (70) 

variance  of  summer  output  z,  (74) 

f itt ing  procedure ,  ( B- 3 ) 

fitting  procedure,  (B-9) 

characteristic  function  of  z,  (G-2) 


j  *  * 

VllX 


TR  10539 


DETECTION  PERFORMANCE  OF  A  MODIFIED  GENERALIZED  LIKELIHOOD 
RATIO  PROCESSOR  FOR  RAND(M  SIGNALS  OF  UNKNOWN  LOCATION 

INTRODUCTION 

The  detection  of  weak  signals  of  unknown  strength  and 
location,  in  the  presence  of  noise,  is  a  coomon  requirement  in 
many  applications.  One  method  of  trying  to  restore  information 
about  the  unknown  parameters  of  the  signal  is  to  estimate  them 
from  the  received  data  and  to  substitute  these  estimates  into  the 
likelihood  ratio.  When  this  is  done,  the  resulting  generalized 
likelihood  ratio  is  compared  with  a  threshold  for  a  decision 
about  signal  presence  or  absence. 

Here,  we  will  consider  a  modification  of  the  generalized 
likelihood  ratio  procedure,  to  be  called  the  MGLR  processor, 
where  only  those  signal  strength  estimates  above  a  minioium 
acceptable  level  are  utilized  [1].  This  results  in  a  processor 
with  nonlinearities  in  each  channel  which  have  breakpoints  in 
their  characteristics.  This,  in  turn,  leads  to  chazmel  outputs 
and  a  processor  output  which  are  distinctly  non-Gaussian  in  many 
cases,  thereby  making  the  analysis  of  performance  a  rather 
difficult  task.  In  fact,  use  of  the  Gaussian  approximation  and 
the  central  limit  theorem  can  yield  very  misleading  values  for 
the  threshold  settings  and  signal-to-noise  ratios  that  are 
required  in  order  to  achieve  specified  false  alarm  and  detection 
probabilities . 
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The  specific  problem  of  interest  is  as  follows.  A  total  - 
number,  N,  of  search  bins  contains  independent  noises  of  uniform 
power  distribution;  this  can  be  accomplished  by  a  normalization 
procedure,  ffhen  signal  is  present,  it  occupies  a  subset  of  M  of 
these  bins,  whose  locations  are  unknown.  Each  of  the  N  channel 
outputs  are  subjected  to  a  common  nonlinearity  and  then  summed; 
this  sum  is  then  compared  with  a  threshold.  We  will  derive  the 
false  alarm  and  detection  probabilities  of  this  M6LR  processor 
and  verify  the  results  by  simulation  over  a  wide  variety  of 
parameter  values. 

Some  related  results  on  a  processor  with  a  dead  zone  are 
available  in  [2];  however,  the  work  here  is  a  generalization  in 
several  respects.  The  number  of  signal  bins,  M,  can  be  less  than 
the  total  number  of  search  bins,  M,  and  the  approximating 
nonlinearity  is  more  general,  being  characterized  by  two 
parameters  rather  than  one.  In  fact,  in  the  following  analysis 
for  the  detection  probability,  we  will  utilize  a  different  pair 
of  parameter  values  for  the  signal-plus-noise  bins,  compared  with 
the  pair  for  the  noise-only  bins.  The  fact  that  M  is  less  than  N 
significantly  complicates  the  detection  probability  analysis 
relative  to  that  in  [ 2 ] . 
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PROBLEM  DEFINITION 

A  search  space  consisting  of  N  bins  contains  a  noise 
background  which  has  a  uniform  power  distribution  over  all  the 
bins;  this  situation  is  achieved  in  practice  by  the  use  of 
normalization,  which  could  be  based  upon  neighboring  bins  or  past 
behavior.  In  addition,  when  signal  is  present,  N  of  these  bins 
contain  random  signal  components  of  unkno«m  locations  and 
strengths.  We  will  presume  that  the  number  of  occupied  bins,  M, 
is  also  unknown;  this  will  lead  to  a  processor  that  must  stake 
some  asstmption  about  M  in  order  to  achieve  decent  levels  of 
performance  in  terms  of  the  false  alarm  and  detection 
probabilities.  Other  approaches,  which  estimate  M,  or  which  can 
function  without  knowledge  of  M,  are  the  subject  of  a  future 
report  [ 3 ] . 

Under  the  noise-only  hypothesis  Hq,  the  N  bin  outputs, 
are  statistically  independent  and  identically-distributed  random 
variables  with  exponential  probability  density  function 

qo(Un)  ■  exp(-Ujj)  for  Ojj  >  0  ^or  1  ^  n  ^  N  .  (1) 

For  example,  these  random  variables,  {z^},  could  be  the  envelope- 
squared  outputs  of  a  bank  of  disjoint  narrowband  filters  (or  FFT 
outputs)  for  a  Gaussian  noise  input  with  a  flat  spectrum  over  the 
total  search  band.  The  assumption  of  a  unit  average  power  level, 
E{Xn}  *  ^  *  J  du  qQ(u)  u  -  1  in  (1),  is  a  matter  of  convenience 
and  is  considered  to  have  been  achieved  in  practice  through  the 
use  of  normalization. 
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Under  the  signal>plus-noi8e  hypothesis  N  of  the  bins  also 
contain  signal,  with  average  power  ^  in  the  n-th  bin.  These 
constants,  unknovm.  The  governing  probability  density 

function  for  the  bin  output  random  variables  {z^}  in  this  case  is 

QlCUn)  -  S„  <»P(-S„  %)  tor  u  >  0  ,  a„  - 

for  n  corresponding  to  the  N  occupied  signal  bins;  otherwise, 

density  (1)  holds  true.  Since  the  actual  n-th  signal  power  level 

is  unknovm,  it  must  be  estimated  from  the  available  data 

measurements  or  observations,  of  the  N  bin  outputs. 

The  derivation  of  the  MGLR  (suxlified  generalized  likelihood 

ratio)  processor  [1]  for  this  environment  is  accomplished  in 

appendix  A.  In  order  to  introduce  it,  we  first  define  as  the 

minimum  acceptable  signal  power  level  that  will  be  tolerated  as 

an  estimate  of  the  signal  strength.  Also,  we  define  breakpoint 

■  1  S^,  and  nonlinear  function 

o  o 

fi  -  1  -  ln(x)  for  I  i  X  1 

9(*)  •  {  •  (3) 

[  0  otherwise  } 

We  require  level  1  0,  thereby  making  k  1.  Then,  the  MGLR 
test  for  deciding  on  presence  or  absence  of  a  sigxial  in  obser¬ 
vation,  (x^),  depends  on  decision  variable  z,  which  is  given  by 

We  are  interested  in  determining  the  false  alarm  and  detection 
probabilities  of  this  MGLR  test,  where  v  is  the  threshold. 
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DERIVATION  OF  FALSE  ALARM  PROBABILITY 


APPROACH 

An  exact  derivation  of  the  false  alara  and  detection 
probabilities  of  the  precise  NGUl  processor  in  (3)  and  (4)  %K>uld 
be  extreiaely  difficult  and  tia«*consusiing.  Instead,  we  will 
determine  the  performance  by  cos^ining  two  different 
approximation  techniques.  In  the  first  technique,  the  breakpoint 
nonlinearity,  g(x),  in  (3)  is  replaced  by  the  best  fitting 
piecewise«l inear  device,  which  is  then  analyzed  exactly.  On  the 
other  hand,  in  the  second  technique,  the  exact  device,  g(x),  is 
retained,  but,  instead,  the  probability  analysis  is  approximate, 
relying  on  a  fourth-order  fit  to  the  exact  first  four  moments  or 
cvuBulants  of  decision  variable  s.  This  two-pronged  approach  has 
been  found  adequate  to  cover  the  range  of  values  of  N,  M,  x^,  and 
signal-to-noise  ratio  that  are  of  major  interest.  Simulations 
verify  the  accuracy  and  viability  of  this  approach. 


PIECEWISE-LINEAR  APPROXIMATION 

The  specified  nonlinearity  g(x)  of  interest  is  given  by  (3). 
However,  the  approximating  nonlinearity  that  we  adopt  in  its 
stead  is  the  piecewise-linear  function 

{c  +  b(x  -  x^)  for  X  1  x_l 

°  ,  (5) 

0  for  X  <  XqJ 


5 


TR  10539 


where  constants  c  and  b  depend  on  breakpoint  nonlinearity 
g(x),  and  the  probability  density  functions  of  input  randcm 

variables  as  given  in  (1)  and  (2)  for  k  •  0  and  1/ 

respectively.  For  example,  we  could  choose  the  tangent  line  at 
breakpoint  x^  (independent  of  g|^(u));  then,  we  %rauld  have 

c  -  gCx^)  -  Xq  -  1  -  InCx^)  ,  b  -  g'(Xjj)  -  1  -  .  (6) 

Other  choices  for  c  and  b,  which  minimize  a  mean  square  error 
criterion,  are  derived  and  tabulated  in  appendix  B.  Additional 
vievrpoints  and  interpretations  are  presented  there  and  should  be 
examined  for  a  complete  explanation  of  the  approach. 

The  approximating-nonlinearity  output  is  no  longer  the  random 
variable  y^^  ■  introduced  in  (4),  but  is  instead  the  random 

var j  able 

hn  *  for  1  S  n  S  N  .  (7) 

The  approximation  to  the  original  summer  output  z  in  (4)  is 


C 


N 


1 


(8) 


We  will  evaluate  the  exceedance  distribu  fon  functions  of  this 
new  random  variable,  exactly  for  both  hypotheses  Hq  and  H^. 
Hopefully,  the  piecewise-linear  fit  utilized  in  (5)  will  be 
adequate  to  accurately  determine  false  alarm  probability  in 
the  10~^  range  and  detection  probability  in  the  .99  range. 
Simulation  results  will  bear  out  this  approach. 
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EXCEEDANCE  PROBABILITY  FOR  NOISE  ALONE 

For  noise  alone,  the  probability  density  function  and 
exceedance  distribution  function  of  input  randoa  variables  (x^) 
to  device  h(x)  are,  from  (1),  res{>ectively, 

QqCu)  «  exp(-u)  for  u  >  0  , 


Qo(v|  ■  J  du  q^Cu)  -  exp(>v)  for  v  >  0  .  (9) 

V 

Parameters  c  and  b  of  the  approximating  nonlinearity  h(x)  are 

positive  constants;  in  particular,  in  the  case  Bq  of  noise-alone 

inputs  {x^},  they  will  be  denoted  by  c^  and  b^,  respectively. 

This  reflects  and  accents  the  fact  that,  during  analysis .  the 

approximating  nonlinearity  h(x)  should  be  chosen  differently  for 

the  noise-only  bins  versus  the  signal-plus-noise  bins  in  case  B^^. 

Ifv>c^>0,  the  exceedance  distribution  function  of  device 
o 

output  -  h(Xjj)  in  (7)  is  (for  b^  >  0) 


Ej^(v)  -  Pr(li„  >  V)  -  Pr(Cj,  +  h^ix^  -  x^,)  >  v)  - 


(10) 


k 


V  -  c. 


Pr  x_  >  x^  +  - r- 

o  b. 


V  -  c, 

I 

o 


where  we  used  ( 5 )  and  ( 9 ) ,  and  where 


B  -  exp(-XQ)  ,  Bj  -  1  -  B  . 


(11) 


The  probability  density  function  of  device  output  follows  as 
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E  *  E"  “Pf'  “  b  °)  *  ®o> 


for  all  u,  where  U  is  the  unit  step  function.  The  characteristic 
function  of  is  therefore  the  Fourier  transfora 


-  Bj  ♦  B 


exp(iec^) 


for  all  I  . 


(13) 


The  approx ination  C  (to  the  original  sussoer  output  s)  is 
given  by  (8)  as  a  sum  of  M  independent  terns  Therefore, 

its  characteristic  function  is  closed  fora 


„  exp(UCo) 

-  *^(0  -  +  B 


-B».^(«)BrB»^ 


exp(ilnCp) 


(1  -  iEb„)" 


(1«) 


At  this  point,  it  is  convenient  to  define  several  sets  of 
functions  which  will  help  to  ainiaize  notation  in  the  sequel. 

The  first  is  the  set  of  nomalized  probability  density  functions 


fu””^  exp(-ul 


Pn(u) 


(n-l)l 


for  u  >  0 


0  for  u  i  0 


for  n  *  1,2,... 


(15) 


with  their  corresponding  characteristic  functions 


f„(5)  »  f  du  exp(i(u)  P„(u)  -  - - - -  for  n  •  1,2,...  (16) 

n  J  n  j  n 

wCD 
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Also*  define  the  corresponding  exceedance  distribution  functions 

m 

B„(v)  -  J  dn  p„(a)  - 

V 

Finally,  define  the  corresponding  cusmlative  distribution 
functions 

Cn(v)  •  1  -  E,j(v)  for  all  v  ,  for  n*l,2,...  (18) 

Then,  the  probability  density  function  of  C,,  corresponding  to 
characteristic  function  (14),  is 

p^(u)  -  bJ  »(u)  +  ^  (J)  B^*”  rf*  ^  Pn(“  bp*””) 

where  we  used  (15)  and  (16).  The  corresponding  exceedance 
distribution  function  of  C  iSr  froa  (17), 

Ej.(v)  -  J  da  p^(u)  -  g  (J)  bJ-"  b"  Bnf  bp)  V  >  0.  (20) 

Finally,  the  cumulative  distribution  function  of  ^  is,  from 
(18)  and  (19), 

Cj;(v)  -  B?  +  ^  (3  Bf"  b"  for  V  >  0  .  (21) 

It  is  very  important  to  observe  in  (21)  that  as  soon  as 
summation  index  n  exceeds  v/c^,  this  sum  can  be  terminated, 
because  cumulative  distribution  function  in  (18)  is  zero  for 


n-1  k 

exp(>v)  J  for  v  >  0 
1  for  V  5  0 


for  n  bb  1,2,... 

(17) 
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negative  argunents.  This  is  a  very  advantageous  property, 
especially  for  large  values  of  breakpoint  because  paraaeter 
c^  is  large  then.  That  is,  the  upper  limit  N  in  (21)  can  be 
replaced  by  the  integer  N'  ■  lllT(v/c^),  which  can  be  considerably 
less  than  N. 

This  result  in  (21)  for  cumulative  distribution  function 
C^(v)  pertains  to  hypothesis  Hq,  noise-only  present  in  all  N 
search  bins.  Parameters  c^  and  b^  used  in  device  h(z)  in  (5)  are 
chosen  for  a  good  fit  to  g(z)  in  (3)  for  z  >  z^  for  the  case  of 
noise-only  inputs.  See  apx>endiz  B  for  more  details  on  the 
fitting  procedure  and  numerical  results. 
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DERIVATION  OF  DETECTION  PROBABILITY 

STATISTICS  OP  BIN  OUTPUTS  WITH  SIGNAL 

Por  the  signal -present  bins,  the  probability  density  function 
of  inputs  (z^)  is  given  by  (2).  The  corresponding  exceedance 
distribution  function  for  these  particular  bins  is 

Qj(v)  »  ezp(-av)  for  v  >  0  .  (22) 

Here,  for  the  analysis  of  perfomance  of  susoser  output  C,  in  (8), 
we  presume  that  the  signal  power  levels  S^  in  all  M  occupied  bins 
are  equal,  with  value  S;  then  a^*a*  (1  ‘tS)"^  for  these  N 
bins,  while  a  «  1  for  the  remaining  N  -  N  noise-only  bins. 

For  the  signal-plus-noise  bins,  we  allow  the  piecewise-linear 
approximation  h(x)  in  (5)  (to  the  given  nonlinearity  g(x))  to 
have  different  values  for  the  parameters,  c  and  b,  which  depend 
on  S,  x^,  and  g(z).  The  values  c^  and  b^  for  the  N  -  N  noise- 
only  bins  are  kept  the  same  as  in  the  false  alarm  analysis  above. 

Ifv>c>0,  the  exceedance  distribution  function  of  a 
signal-present  bin  output  1)^^  is  (for  b  >  0) 

E^(v)  -  Pr(l|n  >  V)  -  Pr(c  +  b(Xjj  -  x^)  >  v)  -  (23) 

■  >  *0  +  ■  »  “P(- 

where  we  used  (22)  and  defined 

A  -  exp(-axQ)  ,  Aj^  ■  1  -  A  .  (24) 
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The  probability  density  function  of  device  output  is 

p^(u)  “  Aj  6(u)  +  a  5  e*p[“  j  U(u  -  c)  for  all  u  ,  (25) 

while  the  corresponding  characteristic  function  is 

f^(t)  -  Aj  +  A  '  •  <**> 

These  functions  in  (25)  and  (26)  pertain  only  to  the  signal-plus- 
noise  bins. 


STATISTICS  OF  SUMMER  OUTPUT  C, 

Since  summer  output  C  is  composed  of  M  (^1)  signal  bins  and 
K  «  M  -  M  noise-only  bins,  its  characteristic  function  is 


where  we  used  (26)  and  (13).  If  S  «  0,  then  a  «  1,  A  •  B, 

Aj  ■  Bj,  c  «  Cq,  b  ■  b^,  and  fj.(^)  reduces  to  the  earlier  result 
in  (14)  for  noise  only. 

For  M2  1,  we  use  the  binomial  expansion  on  (27)  to  get 


t.(E)  ■  (a,  +  A  ■  “o  C  “n  _  ,  (28) 

*  11  1  -  ®  W  "  (1  -  iEb/a)" 


where 


D,  -  0,(H,aXg)  a  g)  A^-"  a"  for  0  S  »  S  M 


(29) 
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For  K  2  1,  that  is,  M  <  N,  (the  case  of  M  •>  N  is  nuch  siapler), 
the  noise-only  characteristic  function  tern  is 


«k(^) 


f  exp(iCc  )»K  K  ezp(i 


ezpCiUc^) 


(1  -  iCb^)- 


.  (30) 


where 


■  (k)  B*'  for  0  i  k  i  K  .  (31) 

Therefore,  the  characteristic  function  of  summer  output  C#  for 
signal  present  in  M  out  of  N  bins,  is 


expdCkc^) 


0  ^  *  (1  -  itb  )*  °  m-1  ®  (1  -  Ub/a)“ 


M  K  exp( i5 (mc+kc_ ) ) 

+  HZ  rz  Un.  - - — 2 - j- . 

m-l  k-1  ®  (1  -  i5b/a)“  (1  -  Ub_)* 


(32) 


At  this  point,  we  utilize  the  results  in  appendix  C  (by 
identifying  a  ->  b/a,  fi  b^,  z  -►  -iC  in  (C-9))  to  obtain  the 
partial  fraction  expansion  of  the  last  term  in  (32): 


5=i  ^  ^ 


(1  -  i£b/a)“  (1  -  iEb-)*^  (1  -  iEb/a) 


D’  ^  1 . .  nm- 

W  (1  -  i 


(1  > 


(33) 


Substitution  in  (32)  then  yields  the  characteristic  function  of 
summer  output  C  in  the  form 
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K  exp(i5kc  )  M  exoritmc) 

C  0  0  0  ^  (1  -  ilh\^  °  bSTl  (1  -  iCb/a)“ 


M  K 


+  C  ^  %  \  exp(ie(«c+kcQ))  X 
(1  -  itb/a)®”"  (1  -  ilb 


(34) 


We  now  use  the  following  connection  between  a  corresponding 
characteristic  function,  its  probability  density  function,  and 
its  cumulative  distribution  function;  referring  to  (15)  -  (18), 
we  have,  for  k  integer. 


exp(iM) 
(1  -  UtV 


?  Pk(H^)  torkil  .  (35) 


The  cumulative  distribution  function  of  C  corresponding  to 
characteristic  function  (34)  then  follows  immediately  in  the  form 


1  Co(v)  + 

K  #v-kc  . 

^0  ’^k  ^k(  b  , 

k«l  o 

ZZ  <'kn<‘>''a'‘>o> 

>1 

ln“0 

k-1 

,v-mC’ 

+  C  >»^(b/a,b^) 

,v-mc-kc. 


'))■ 


(36) 


where  Cq(v)  «  0  for  v  <  0,  and  is  1  otherwise.  One  of  the 
attractive  features  of  this  expression  is  that  when  mc+kc^ 
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exceeds  threshold  v,  normalized  cumulative  distribution  functions 
^m>n  ^k-n  zero,  allowing  the  summations  on  m  and  k  to  be 
terminated. 

Coefficients  and  in  (34)  and  (36)  are  immediately 
available  from  appendix  C  by  identifying  a  b/a,  fi  -*  b^  in  (C-1) 
and  (C-2)  to  get 


^kn 


mn 


(I'-o) 


b*'  (-Sb„)"  (lc)„ 


(b-ab  ) 


k-fn  nl 


for  k21,n20. 


(37) 


(ab_)“  (-b)“  (m), 


(ab  -b) 


m-t'n  nl 


for  m21,nil0.  (38) 


These  quantities  can  be  rapidly  evaluated  by  recursions;  dropping 
the  argiments  for  the  time  being,  we  have 

"ko  -  (b^b;)''  *or  k  k  1  ,  (39) 

"kn  -  »'k,n-l  ^  ^  for  n  k  1  ,  k  k  1  ,  (40) 


along  with 


(41) 


m-l+n 


mn 


m,n-l  b-ab 


for  n  2  1  ,  m  2  1 


(42) 


The  desired  result  for  the  cumulative  distribution  function 
of  summer  output  C  is  given  by  (36).  The  detection  probability 
for  this  processor  is  1  -  C^(v),  for  threshold  v. 
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FOURTH-ORDER  MOMENT  APPROACH 

When  breakpoint  in  nonlinear  device  9(x)  in  (3)  is  large, 
the  techniques  in  the  previous  two  sections  %#ork  well  enough  to 
yield  useful  results  for  both  the  false  alarm  and  detection 
probabilities.  However,  as  x^  decreases,  the  series  expansions 
require  progressively  more  terms,  and  the  execution  time  becosies 
excessive  for  large  values  of  M,  like  1024.  In  addition,  the 
alternating  series  in  (36)  for  the  cumulative  distribution 
function  loses  all  its  significance,  eventually  giving  negative 
values  as  well  as  values  larger  than  1. 

Another  approach  that  has  worked  well  in  the  past  is  to  use 
a  large  nxunber  of  moments  or  cusiulants  of  the  decision  variable 
of  interest,  if  they  can  be  calculated  accurately  to  high  order 
[4].  Additional  elaboration  on  this  approach  is  given  in 
appendix  D.  Here,  however,  we  will  limit  consideration  to  the 
first  four  moments  of  the  exact  decision  variable  s  in  (4)  and 
fit  a  fourth-order  candidate  characteristic  function;  this 
procedure  has  been  used  successfully  in  a  previous  detectability 
study  [5].  We  label  this  procedure  FIT-4  here. 

Let  the  k-th  cumulant  of  random  variable  ^ 

denoted  by  Xy(l^/3c^,S) ,  where  it  is  presumed  that  the  signal 
powers  in  the  occupied  bins  are  all  equal  to  S.  (The  dependence 
of  the  cumulants  on  breakpoint  x^  is  indicated  explicitly. )  For 
the  noise-only  bins,  the  pertinent  cumulant  is  . 

When  signal  is  present  in  M  bins  out  of  the  total  of  N  search 
bins,  the  k-th  cumulant  of  summer  output  z  in  (4)  is  given  by 
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XjCJt)  •  M  Xy(k,XQ,S)  +  (N  -  M)  Xy(k,Xjj,0)  for  k  i  1  .  (43) 

Xn  order  to  find  the  cumulants  IXy(l^f<o#S) ),  we  first  nunerically 
evaluate  the  corresponding  moments  (Pylk^x^^S))  exactly  by  means 
of  the  integral  expression 

Py(lt»aCo'-^  “  ^  “  J  9l(«)  9(u)*^  - 

m 

«  J  du  a  exp(-au)  (u  -  1  -  ln(u))^  for  k  l  1  ,  (44) 

*o 

along  with  Py(0,x^«S)  •  1.  In  particular,  we  evaluate  (44)  for 
k  •  1,2, 3, 4.  The  desired  cumulants  of  y^  then  follow  according 
to 

Xy(l»*o'-^  ■  ' 

Xy(2,XQ,S)  -  #iy(2,XQ,S)  -  /iy(l#*o/S)  , 

Xy(3)  -  #iy(3)  -  3  /iy(2)  #iy(l)  +  2  ,  (45) 

Xy(4)  -  /iy(4)  -  4  /iy(3)  py(l)  -  3  +  12  fiy(2)  -  6/iJ(l) 

where  we  have  adopted  an  abbreviated  notation  in  the  last  two 
lines.  Tables  of  the  first  four  cumulants  (XyCl^r^o'S) )  are  given 
in  appendix  E. 

The  candidate  characteristic  function,  FIT>4,  that  will  be 
fit  to  the  first  four  cumulants  {X2(^)}  (43)  is  [5;  page  59] 
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£4(^)  *  SSPlASb*  ~  -  expfiCb'  -  x}c* n  -  M'  ln(l  -  i^w')] 

*  (1  -  i^w*)"  ^  * 

(46) 

which  is  a  product  of  Gaussian  and  chi-squared  characteristic 
functions;  this  guarantees  that  f^CC)  is  a  legal  characteristic 
function,  although  its  corresponding  probability  density  function 
can  be  nonzero  for  negative  arguments.  For  notational  shorthand, 
let 

*1  -  X*(l)  ,  X2  -  X*(2)  ,  X3  -  X*(3)/2  ,  X4  -  X,(4)/6  .  (47) 


Then,  in  order  for  characteristic  function  (46)  to  have  cusiulants 
(X2(X))  for  k  •  1,2, 3, 4,  we  must  take  its  constants  as 


b' 


(48) 


2 

If  X2  <  X3/X4,  then  c*  <0;  in  this  case,  an  alternative  is 
to  set  c'  «  0,  getting  instead  the  constants 


b' 


(49) 


This  is,  in  fact,  the  third-order  fit  utilized  in  [5;  page  43]. 
However,  even  if  c'  <  0  in  (48),  candidate  characteristic 
function  (46)  can  still  be  useful;  namely,  the  magnitude  |f4(C)i 
has  a  minimum  at 


(50) 
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with  mininum  value 

l*4(t,)l  -  e.p(f- (f^]  .  (51) 

If  ^4(^1)  I  <«  <4(0)  *  1*  namely  rapid  decay  of  f4(C)#  %re  can 
still  use  approximation  f^(C)  in  (46)  to  get  useful  results  for 
the  corresponding  exceedance  distribution  function,  by  utilising 
only  those  values  K|  < 
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ENERGY  DETECTOR  PERFORMANCE 

This  section  is  not  concerned  with  perfonaance  of  the  NGLR 
processor  or  its  nonlinearity  9(z)  given  by  (3).  Rather,  we 
derive  the  exact  performance  of  a  processor  esqploying  a 
completely  linear  summation  of  the  same  N  search  bin  outputs 
(z^),  looking  for  a  signal  in  M  bins  of  unknown  location.  This 
system  will  serve  as  a  worthwhile  basis  for  comparison  with  the 
MGLR  processor  and  will  indicate  quantitatively  when  we  should 
prefer  one  processor  over  the  other.  Some  related  results  for 
weighted  energy  detectors  operating  under  missiatches  in  frequency 
and  time  locations  are  available  in  [6]. 

The  energy  detector  sums  over  all  N  bin  outputs  according  to 
the  linear  sum 


N 


It  is  called  an  energy  detector  because  bin  outputs,  nre 

themselves  envelope-squared  quantities.  By  reference  to  (1)  and 
(2),  it  is  immediately  seen  that  the  characteristic  function  of  w 
under  hypothesis  is 


(1  -  it/a)"  (I  -  U)"***  ' 


a  • 


1 

i  +  s  ' 


(53) 


where  we  presume  equal  signal  power,  S,  in  all  N  occupied  bins. 
At  this  point,  we  could  directly  use  the  partial  fraction 
expansion  presented  in  appendix  C;  however,  we  would  again 
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encounter  an  alternating  series  and  loss  of  significance. 

Instead f  we  have  resorted  to  the  accurate  nusierical  FFT  procedure 
presented  in  [7]  in  order  to  determine  the  receiver  operating 
characteristics  exactly. 

Sone  special  cases  of  (53)  are  %#orth  noting.  For  signal 
power  S  •  0,  the  characteristic  function  reduces  to 


1 

II  » 

(1  -  U)** 


(54) 


for  which  the  false  alarm  probability  is 


N-1  n 

-  Pr(w  >  v|Hq)  -  exp(-v)  ^  ^  for  v  >  0  . 


(55) 


Another  special  case  occurs  for  M  ■  N,  namely  signal  in  every 
bin,  when  present.  The  energy  detector  is  then  the  optisnim 
(lilcelihood  ratio)  processor.  The  characteristic  function  in 
(53)  reduces  to 


VO 


1 

(1  -  U/a)**  ' 


a 


1 

TT~s  ' 


(56) 


for  which  the  detection  probability  is 


P^  -  Pr(w  >  v|Hj) 


exp(-av) 


for  V  >  0  , 


(57) 


while  the  false  alarm  probability  is  (55)  again. 
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6RAPBICAL  RESULTS 
PIBCBWISB-LINBAR  FIT  h(x) 

An  example  of  the  ability  of  piecewise-1 inear  fit  h(x)  in 
(5),  to  approxinate  the  actual  nonlinearity  g(x)  in  (2),  is  given 
in  figure  1  for  signal  absent,  S  «  0,  and  for  breakpoint  values 
x^  *  1,3,5,7,9,11.  As  expected,  the  approxiaations  for  the 
larger  x^  values  are  on  top  of  their  corresponding  portions  of 
the  dashed  curve,  which  is  x  -  1  -  ln(x).  However,  for  the 
smaller  values  of  breakpoint  x^,  like  1,  the  best  linear  fit  cuts 
across  the  dashed  curve,  and  for  large  x  values,  considerably 
underestimates  the  desired  values.  This  is  unavoidable  because 
of  the  considerable  curvature  of  g(x)  near  x  *  1  when  x^  ■  1. 

The  fit  is  significantly  improved  for  •  3  and  larger  values. 

The  example  in  figure  2  holds  the  breakpoint  fixed  at  value 
x^  ■  5  while  the  signal  power,  S,  is  varied  over  the  values  0  dB, 
5  dB,  10  dB.  The  fit  for  S  >  5  dB  lies  on  top  of  the  dashed 
curve,  while  the  fit  for  S  •  10  dB  overestisiates  g(x)  in  the 
plotted  region  displayed.  However,  for  argusient  values  x  larger 
than  18,  the  two  curves  will  eventually  cross. 

These  results  in  figures  1  and  2  indicate  that  we  can  expect 
the  piecewise-linear  fitting  procedure  to  yield  useful 
performance  prediction  results  for  the  larger  values  of 
breakpoint  x^.  On  the  other  hand,  they  also  indicate  that  some 
problems  can  be  expected  for  the  smaller  values  of  x^,  and  that 
an  alternative  procedure  will  be  necessary  for  those  cases. 
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FALSE  ALARM  PROBABILITY 

The  false  alarm  probability  for  the  M6LR  processor  was 
evaluated  numerically  from  the  cusnilative  distribution  function 
in  (21),  namely  as  1  -  C^(v),  for  various  values  of  search  size 
N  and  breakpoint  x^,  versus  threshold  v.  These  results  are 
presented  in  figures  3-7  for  M  «  1024,  512,  256,  128,  64, 
respectively.  For  example,  in  figure  3  for  N  >  1024,  breakpoint 
is  varied  from  1  to  11  in  increments  of  1.  Superposed  on 
these  smooth  curves  in  figure  3  are  a  set  of  jagged  curves 
corresponding  to  a  simulation  utilizing  700,000  trials  of 
processor  output  z  using  the  exact  nonlinear  device  g(x)  in  (4). 
This  simulation  verifies  the  accuracy  of  (21)  over  the  coiq>lete 
range  of  values  in  figure  3,  including  the  low  threshold  values, 
V,  in  the  range  (0,10),  where  the  exceedance  distribution 
function  has  2d)rupt  changes  in  its  character. 

Similar  confirmation  is  indicated  in  the  remaining  figures 
4-7,  where  simulations  of  various  sizes  (at  least  1,000,000 
trials)  have  been  employed  for  comparison.  The  only  case  of 
significant  departure  between  theory  and  simulation  is  in  figure 
7  for  small  N  (64),  small  x^  (<  3),  and  large  v  (>  20).  The 
systematic  deviation  between  the  two  curves  for  x^  >  2  (at  small 
Pf  values)  illustrates  that  the  piecewise-1 inear  fit  described  in 
appendix  B  is  unable  to  accurately  represent  device  g(x)  in  this 
particular  region. 
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The  particular  results  in  figures  3-7  were  obtained  by 
means  of  FIT-2,  which  is  derived  and  tabulated  in  appendix  B. 
That  is,  both  constants  c  and  b  in  linear  fit  (5)  were  optimized 
for  minimum  mean  squared  error*  Other  fits  fared  poorer  than 
FIT-2  in  the  problem  region  of  figure  7.  For  example,  FIT-1  was 
never  better  than  FIT-2  and  was  sometimes  poorer  in  predicting 
values.  Attempts  to  utilize  a  Gaussian  approximation  for 
decision  variable  z  were  much  poorer,  often  leading  to  gross 
underestimates  of  the  actual  false  alarm  probabilities.  Finally, 
when  FIT-4  in  (46)  was  tried,  it  failed  because  the  fitting 
constants  b'  and  c'  in  (48)  were  negative,  and  remedy  (51)  was 
not  small  enough  to  warrant  its  use. 

In  order  to  rectify  the  shortcoming  in  the  bottom-right 
portion  of  figure  7,  it  would  be  necessary  to  employ  a  two-piece 
linear  fit  to  g(x).  Nevertheless,  as  it  is,  the  combination  of 
the  two  types  of  results  available  in  figures  3-7  enable 
accurate  selection  of  threshold  values,  v,  in  order  to  realize 
specified  false  alarm  probabilities  >  lE-6  for  various 
breakpoint  values  x^  and  search  sizes  N. 
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Figure  3.  versus  v  for  N  ■  1024,  FIT-2 
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V 

Figure  6.  versus  v  for  N  -  128,  FIT-2 
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DETECTION  PROBABILITY 

A  number  of  receiver  operating  characteristics  (ROCs)  for  the 
MGLR  processor  are  presented  in  figures  8-43.  They  have  all 
been  evaluated  for  search  size  N  *  1024.  The  values  of  the 
number,  N,  of  bins  occupied  by  signal  and  the  breakpoint,  z^, 
cover  all  36  combinations  f  the  following  numbers: 

M  -  8,  16,  32,  64,  128,  256, 

Zq  -  1,  3,  5,  7,  9,  11.  (58) 

A  detailed  ezplanation  of  the  ROC,  receiver  operating 
characteristic,  in  figure  8  for  N  >  1024,  M  •  8,  z^  •  1,  axul 
FIT-4  follows.  The  solid  smooth  curves  were  evaluated  by  means 
of  the  candidate  characteristic  function  in  (46),  for  both  the 
false  alarm  probability  as  well  as  the  detection  probability. 

The  jagged  curve,  labeled  A,  is  a  simulation  result,  using  10,000 
trials  with  the  ezact  nonlinear  device,  g(z),  with  no 
approzimations ,  for  S  ■■  9  dB.  The  simulation  result  is  seen  to 
overlay  and  verify  the  theoretical  curve  within  random  errors 
inherent  with  a  limited  number  of  trials.  The  dashed  curve, 
labeled  B,  is  what  the  Gaussian  approzimation  predicts  for  the 
performance  capability,  again  for  S  ■  9  dB.  For  this  particular 
example  of  small  z^,  namely  1,  the  Gaussian  result  is  optimistic 
by  only  .6  dB  at  the  left  edge  of  the  curve. 

When  breakpoint  z^  is  increased,  the  behavior  of  the  receiver 
operating  characteristics  for  the  MGLR  processor  changes 
significantly,  and  it  is  necessary  to  switch  to  procedure  FIT-2, 
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which  utilizes  piecewise- linear  function  h(z)  in  (5);  this  switch 
occurs  in  figure  10  for  2  5  when  N  ■  8.  For  still  larger 
breakpoint  values,  such  as  •  11  in  figure  13,  there  are 
pronounced  bumps  in  the  characteristics,  in  addition  to  an 
unachievable  region.  That  is,  the  only  value  of  false  alarm 
probability,  P^,  larger  than 

**fm  ■  ^  ”  exp(-XQ))**  -  .0170  (59) 

that  can  be  realized  is  1,  due  to  the  dead  zone  character  of 
nonlinearity  g(x).  This  is  due  to  the  impulse  of  area  at  the 
origin  of  probability  density  function  P^(u)  in  (12). 

The  simulation  result  A  for  10,000  trials  and  S  ■  8  dB  in 
figure  13  confirms  the  theoretical  approach,  7IT-2,  to 
performance  prediction.  Another  point  to  observe  is  the  poor 
quality  of  the  Gaussian  approximation  B,  which  cuts  across  the 
plot  with  virtually  zero  slope.  Curve  B  bears  no  resemblance  to 
the  actual  characteristics,  although  it  does  cross  the  correct 
curve  just  below  -  .001. 

Most  of  the  theoretical  receiver  operating  characteristics  in 
figures  8-43  have  been  verified  by  simulations,  each  using 
10,000  trials.  In  no  case  is  the  discrepancy  greater  than  .1  dB. 
One  case  was  pursued  further;  namely,  figure  17  used  200,000 
trials  at  S  >  5  dB,  thereby  allowing  a  check  down  to  P^  *  lE-5. 
The  two  results  for  S  »*  5  dB  in  figure  17  are  overlays,  including 
the  blimps  near  P^  equal  .25  and  .07. 
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From  figures  8-43,  it  is  now  possible  to  determine  the 
signal  power  per  bin,  S,  required  for  the  NGLR  processor,  in 
order  to  realize  a  specified  operating  point  in  terms  of  the 
false  alarm  and  detection  probability  pair  denote  the 

pair  10'^,. 5  as  a  low-quality  operating  point,  and  the  pair 
10~^,.9  as  a  high-quality  operating  point. 

In  table  1,  the  required  signal  powers  per  bin,  S,  for  the 
low-quality  point  are  presented  for  all  36  combinations  of  N  and 
Xq  in  (58)  and  for  N  *  1024.  (The  column  headed  ED  is  for  the 
energy  detector,  (52);  its  receiver  operating  characteristics 
will  be  presented  in  a  later  subsection. )  It  is  imsMdiately 
observed  that  for  each  choice  of  M,  there  is  a  best  brealqioint 
value  Zq  to  use  in  the  MGLR  processor,  in  order  to  operate  with 
minimum  signal  power  S.  For  example,  for  M  *  8,  the  best  z^  is 
11,  and  the  minimum  signal  power  is  S  ■  6.3  dB.  This  is  2.5  dB 
better  than  the  signal  power  required  for  x^  •  1,  naisely  8.8  dB. 
It  is  also  4.9  dB  better  than  the  energy  detector,  which  requires 
11.2  dB  for  M  -  8. 

As  M  increases,  the  best  breakpoint  z^  decreases,  eventually 
reaching  its  minimum  value  of  1  for  N  ■  128.  For  N  ■■  256,  the 
MGLR  and  energy  detector  require  the  same  signal  power  level  of 
-4.0  dB.  As  M  tends  closer  to  N,  the  energy  detector  will 
outperform  the  MGLR;  in  fact,  the  energy  detector  is  the  optimum 
processor  for  M  «  N. 

The  corresponding  results  for  the  high-quality  operating 
point  are  given  in  table  2.  Basically  the  same  conclusions 
follow,  although  all  the  signal  powers  are  necessarily  larger 
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than  in  table  1.  The  only  real  difference  is  that  the  best  H6LR 
processor  still  outperforms  the  energy  detector  at  N  •  256, 
requiring  .3  dB  less  signal  power.  This  trend  is  expected  to  be 
reversed  for  M  somewhat  larger  than  256. 


Table  1.  S(dB)  Required  for  N  «  1024,  ■  10*^,  «  .5 


1: 

1 

3 

5 

7 

9 

11 

ED 

8 

8.8 

8.7 

7.8 

7.0 

6.4 

6.3 

11.2 

16 

6.2 

6.1 

5.4 

5.0 

4.8 

5.0 

8.1 

32 

3.7 

3.6 

3.3 

3.2 

3.3 

3.7 

5.0 

64 

1.2 

1.2 

1.3 

1.6 

2.0 

2.6 

2.0 

128 

-1.3 

-1.1 

-0.6 

0.0 

0.7 

1.5 

-1.0 

256 

-4.0 

-3.6 

-2.6 

-1.5 

-0.5 

0.4 

-4.0 

512 

-7.0 

1024 

-10.0 
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Probability  of  False  Alarm 
Figure  ROC  for  N«1024,  M*8,  Xo«3,  FIT--4 
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R:  Simulation,  10000  trials,  dB 

B:  Gaussian  approximation,  dB 
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Probability  of  False  Alarm 
Figure  10.  ROC  for  N"1024,  M“8,  Xo*5,  FIT-2 
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E-6  E-5  E~<4  .  001  .01  .02  .  05  .1  .2  .  3  .  4  .5 

Probability  of  False  Rlarm 
Figure  ll.  ROC  for  N»!024,  M»=8,  Xo=7,  FIT-2 
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Probability  of  False  Alarm 
Figure  12.  ROC  for  N«1024.  M=8,  Xo=3,  FIT-2 
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Probability  of  False  Alarm 
Figure  13.  ROC  for  N-1024,  M-8,  Xo-11,  FIT-2 
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Probability  of  False  Alarm 
Figure  H.  ROC  for  N*1024,  M«1S,  Xo“l,  FIT-4 
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Figure  /S',  ROC  for  N«=1024,  M«16,  Xo^S,  FIT-4 
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Probability  of  False  Alarm 
Figure  16.  ROC  for  N*=1024,  M®1€,  Xo*=5,  FIT-2 
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Figure  H.  ROC  for  N«=1024,  M«16,  Xo«7,  FIT-2 
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Probability  o-F  False  Rlarm 


Figure  18.  ROC  for  N=1024,  M=16,  Xo=9,  FIT-2 
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Probability  of  False  Alarm 
Figure  21.  ROC  for  N-1024,  M«32,  Xo«3,  FIT-4 
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Probability  of  False  Alarm 
Figure  Zl.  ROC  for  N*1024,  M*=32,  Xo“5,  FIT-2 


51 


Probabf  Hty  of 


TR  10539 


Probability  of  False  Rlarm 
Figure  23.  ROC  for  N*1024,  M=*32,  Xo“7,  FIT-2 
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fl:  Simulation,  10000  trials,  dB 
B:  Gaussian  approximation,  5*4  dB 
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Probability  of  False  Alarm 
Figure  24.  ROC  for  N*1024,  M«32,  Xo-9,  FIT-2 


53 


Probability  of  Detection 


TR  10539 


Probability  of  False  fllarm 
Figure  15.  ROC  for  N-1024,  M»32,  Xo-ll,  FIT-2 
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Ppobabflity  of  False  Alarm 
Figure  26.  ROC  for  N»1024,  M=64,  Xo=l,  FIT-4 
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fl:  Simulation,  10000  trials,  §**0  dB 
B:  Gaussian  approximation,  S*>0  dB 
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Probability  of  False  Alarm 
Figure  27.  ROC  for  N**1024,  M*G4,  Xo=3,  FIT-4 
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Probability  of  False  Alarm 
Figure  28,  ROC  for  N«1024,  M*64,  Xo*5.  FIT-4 
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Probability  of  False  Alarm 
Figure  21.  ROC  for  N«1024.  M«64.  Xo»7,  FIT-2 
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R;  Simulation,  10000  trials,  dB 


B:  Gaussian  approximation,  dB 
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Probability  of  False  Rlarm 
Figure  30.  ROC  for  N*1024,  M=G4,  Xo=9,  FIT-2 
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fl;  Simulation,  10000  trials,  dB 
B:  Gaussian  approximation,  S^3  dB 
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Probability  of  False  Rlarm 
Figure  31.  ROC  for  N=1024,  M=64,  Xo=ll,  FIT-2 
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Figure  32.  ROC  for  N=I024,  M=128,  Xo=l,  FIT-4 
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R:  Simulation  10000  trials,  S*“2  dB 
B:  Gaussian  approximation,  S*--2  dB 
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Figure  33.  ROC  for  N=1024,  M=128,  Xo=3,  FIT-4 
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Probability  of  False  Rlarm 
Figure  K  ROC  for  N«1024,  M«128,  Xo«5,  FIT-4 
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Probability  of  False  Rlarm 
Figure  3S:  ROC  for  N»1024.  M«128,  Xo-7,  FIT~2 
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Probability  of  False  Alarm 


Figure  SC.  ROC  for  N«1024,  M«128,  Xo«9,  FIT-2 
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R:  Simulation,  10000  trials,  S«2  dB 
B:  Gaussian  approximation,  §■‘2  dB 
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Figure  37.  ROC  for  N«1C24,  M*=128,  Xo«ll,  FIT-2 
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Figures®.  ROC  for  N*1024,  M=25G,  Xo=l,  FIT-4 
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fl:  Simulation  10000  trials,  S®~5  dB 
B:  Gaussian  approximation,  dB 
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Figure  3^.  ROC  for  N*1024,  M=25G,  Xo=3,  FIT-4 
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fl:  Simulation  10000  trials,  S*“3  dE 
B:  Gaussian  approximation,  5^-3  dB 
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Figured.  ROC  for  N=1024,  M=256,  Xo=5,  FIT-4 
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Figure -tl.  ROC  for  N'=1C24,  M*=25S,  Xo«7,  FIT-2 
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Figure -^2.  ROC  for  N=1024,  M=256,  Xo=9,  FIT-2 
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Figure  43.  ROC  for  N«1024,  M=256,  Xo=ll,  FIT-2 
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VARIATION  OF  PERFORMANCE  WITH  N 

All  of  the  receiver  operating  characteristics  in  figures 
8-43  above  were  computed  for  search  size  N  «  1024.  In  order 
to  ascertain  the  dependence  on  N,  without  the  same  large-scale 
computational  investigation,  we  have  computed  the  receiver 
operating  characteristics  for  N  >  512  and  N  «  256,  but  only 
for  a  few  selected  cases  of  M  and  z^,  namely 


8,  -  11  , 

(60) 

32,  Xjj  -  7  , 

(61) 

256,  Zq  -  1  . 

(62) 

These  particular  pairs  were  selected  because  they  cover  a  wide 
range  of  values  of  M,  and  they  represent  the  best  choices  of 
breakpoint  z^  for  each  particular  N;  see  tables  1  and  2. 

The  receiver  operating  characteristics  for  (60)  are  given  in 
figures  13,  44,  and  45;  those  for  (61)  are  given  in  figures  23, 
46,  and  47;  while  those  for  (62)  are  given  in  figures  38,  48,  and 
49.  Simulations  verify  the  theoretical  analyses  in  every  case. 

The  signal  powers  required  in  order  to  realize  the  two 
operating  points  of  interest  are  given  in  tables  3  and  4, 
respectively.  As  N  changes  by  a  factor  of  2,  the  differences  in 
required  signal  powers  are  0.2  dB  for  M  ■  8  in  (60);  0.5  dB  for 
M  s  32  in  (61);  and  1.1  dB  for  M  «  256  in  (62).  These 
differences  hold  for  both  the  low-quality  and  the  high-quality 
operating  points.  A  very  crude  approzimation  to  these 
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observations  is  H  log(N)  -  %.  These  results  indicate  that  the 
cost  of  having  to  search  a  larger  region  is  not  significant,  at 
least  for  these  sets  of  parameter  values. 


Table  3.  S(dB)  Required  for  ■  10“^,  ■  .5 


M 

*0 

M  -  1024 

N  -  512 

N  •  256 

8 

11 

6.4 

6.2 

6.0 

32 

7 

3.2 

2.7 

2.2 

256 

1 

o 

. 

1 

-5.2 

-6.3 

Table  4.  S(dB)  Required  for  Pj  ■  10*®,  P^  ■  .9 


N 

*o 

N  -  1024 

N  -  512 

N  -  256 

8 

11 

10.6 

10.3 

10.2 

32 

7 

6.1 

5.6 

5.2 

256 

1 

-1.1 

-2.2 

-3.2 
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Probability  of  False  Rlarm 


Figure  "H.  ROC  for  N«512,  M«8,  Xo*ll,  FIT~2 
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Figure  Is.  ROC  for  N«=256,  M«8,  Xo-11,  FIT-2 
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Figure  *46.  ROC  for  N*5!2,  M=32,  Xo*7,  FIT-2 
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>9ure47.  ROC  for  N-256,  M=32,  Xo-7,  FIT-2 
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Figure -48.  ROC  for  N«512,  M«256,  Xo«l,  FIT-4 
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Probability  of  False  R1 arm 
Figure  ROC  for  N-25B,  M*25e,  Xo«l,  FIT-4 
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ENERGY  DETECTOR  CHARACTERISTICS 

The  energy  detector  was  described  in  (52),  and  its 
characteristic  function  was  presented  in  ( 53 ) .  The  false  alara 
probability,  given  by  (55),  is  plotted  in  figure  50  for  search 
sizes  N  ■  64,  128,  256,  512,  1024  for  probability  values  do%m  to 
the  10~^  level.  The  receiver  operating  characteristics  were 
evaluated  exactly,  by  the  methods  in  [7],  and  are  presented  in 
figures  51  -  58  for  N  «  1024  and  N  ■  8,  16,  32,  64,  128,  256, 
512,  1024,  respectively.  The  values  of  required  signal  power  in 
tables  1  and  2  under  the  ED  (energy  detector)  column  were 
obtained  from  figures  51  >  58.  Since  these  receiver  operating 
characteristics  were  obtained  by  an  exact  analysis  of  the  exact 
linear  summation  in  (52),  there  is  no  need  to  resort  to  a 
simulation  or  a  Gaussian  approximation;  hence,  there  are  no 
overlays . 
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Probability  of  False  filarm 
Figure  S’!.  Energy  Detector  ROC  for  N=1024,  M=8 
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Probability  of  False  Rlartn 
Figure  52.  Energy  Detector  ROC  for  N=1024,  M=1S 
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Figure  53.  Energy  Detector  ROC  -For  N=1024,  M*32 
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Probability  of  False  Rlarm 
Figure  Energy  Detector  ROC  for  N*=1024,  M«64 
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Figure  5'5'.  Energy  Detector  ROC  for  N®1024,  M**128 
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Probability  of  False  Rlarm 
Figure  Energy  Detector  ROC  for  N**1024,  M=256 
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Figure  5^.  Energy  Detector  ROC  for  N*»1024,  M«512 
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Probability  of  False  Alarm 
Figure  5R  Energy  Detector  ROC  for  N“1024, 
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EXTENSION  TO  UNEQUAL  SIGNAL  POWERS 

All  the  earlier  results  presusied  that  the  N  signal  bins 
contained  exactly  the  sane  signal  ponrer  S.  In  practice#  the 
values  of  the  signal  powers  in  the  various  bins  will  be 
different.  Instead  of  having  to  determine  receiver  operating 
characteristics  for  each  particular  set  of  unequal  signal 
strengths  (S^l #  we  will  now  determine  effective  values  to  replace 
S  and/or  M#  that  can  be  used  with  the  existing  results. 

Let  the  signal  power  in  the  m-th  signal  bin  be  ^  for 
1  ^  m  1  M.  We  first  consider  an  equivalent  or  effective  signal 
power  ^  to  be  used  in  place  of  the  set  {^).  The  number  of 
signal  bins  is  kept  at  value  M  in  this  particular  approach. 

An  alternative  approach#  where  an  effective  number#  N^#  of 
signal  bins  is  defined  instead#  while  the  signal  power  set 
is  replaced  by  its  peak  value#  is  considered  next.  Finally#  a 
two-parameter  extraction  procedure#  where  both  ^  and  are 
simultaneously  determined#  is  presented. 
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MATCHING  OF  MEAN  OUTPUTS 

For  breakpoint  nonlinearity  9(z)  ■  z  -  1  -  ln(z)  for  z  >  z^, 
define  the  following  mean  function,  with  a  •  1/(1  '**5), 

P(S,Zq)  ■  J  du  a  ezp(-au)  (u  -  1  -  ln(u)]  - 
*o 

-  ezpc-az^)  -  1  + 

A  plot  of  this  function  is  given  in  figure  59.  The  asynptotic 
behavior  for  large  S  is  approzimately  linear,  being  given  by 
[8;  5.1.11] 

/^(S.Zq)  ~  S  -  ln(l  +  S)  +  y  -  *0(5*0  -  la(*o))''S  as  S  -►  •  .  (64) 

Without  loss  of  generality,  assusie  that  the  signal  is  in  bins 
[1,M].  The  mean  of  the  output  of  nonlinearity  g(z),  for  the 
m-th  signal  input,  is  p(^,Zg).  The  corresponding  siean  of  summer 
output  z  is  ezactly 

M 

#^2  -  HZ  +  (N  -  M)  p(0,Zj,)  .  (65) 

m*l 

On  the  other  hand,  if  all  the  sigixal  powers,  {^),  are 
replaced  by  an  effective  signal  power  value,  S^,  over  all  the  M 
signal  bins,  we  obtain  mean  value 

#^2  -  M  #i(S^/*o)  +  (N  -  M)  p(0,Zjj)  .  (66) 

If  we  match  these  two  ezpressions  in  (66)  and  (65)  for  the  mean 
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output,  we  obtain  the  following  equation  that  ^  mist  satisfy: 

<'(Se'*o>  “  H 

Once  the  right-hand  side  of  (67)  has  been  computed  from  the 

specified  M,  (^1*  and  the  required  to  realize  equality  in 

(67)  can  be  found  from  figure  59.  The  use  of  closed  form  (63) 

makes  the  niunerical  summation  in  (67)  a  simple  computer  task, 

once  a  routine  for  exponential  integral  Bj^(z)  is  available. 

If  all  (^)  are  equal  to  S,  then  (67)  yields  ^  *  S,  as 

expected.  However,  in  general,  the  quadratic  behavior  of  mean 

function  near  S  >  0  indicates  that  the  smaller  signal 

powers  are  suppressed  relatively,  the  store  so  for  larger  x^. 

The  example  M  >  4,  «  m  for  1  6  m  1  M,  yields  effective  signal 

powers  ^  «  2.554,  2.557,  2.586,  2.638,  2.700,  2.763  for 

Xq  ■  1,  3,  5,  7,  9,  11,  respectively. 

This  replacement  approach  of  (S^)  by  takes  account  of  the 

specific  nonlinearity  g(x),  the  breakpoint  value  x^,  and  all  of 

the  individual  signal  powers  {^)  for  1  6  m  1  M.  However,  it 

does  require  the  solution  of  (67)  for  argument  ^  which  leads  to 

equality  with  the  computed  right-hand  side. 

Instead  of  keeping  M  fixed  throughout  this  replacement,  we 

could  instead  define  an  effective  nusiber,  M.,  to  be  used  with  the 

e 

replacement  of  the  set  of  signal  powers,  (S||^),  by  its  peak  value 
Sp  «  max{^).  In  this  latter  case,  the  first  term  on  the  right- 
hand  side  of  (65)  becomes  fj(Sp,x^),  When  we  equate  this 
expression  to  the  corresponding  term  in  (65),  we  have  an  explicit 


M 
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result  for  the  effective  number  M_  of  signal  bins  as 


-e; 

m*l 


(68) 


This  result,  in  conjunction  with  expression  (63),  allows  for 
rapid  and  simple  evaluation  of  M^,  once  and  set  {^)  are 
specified.  For  the  same  numerical  example  above,  namely  M  ■  4, 
s  ■<  m  for  1  ^  m  ^  M,  we  find  ^  *  4  and  N.  *  2.409,  2.343, 
2.184,  2.020,  1.876,  1.756  for  x^  -  1,  3,  5,  7,  9,  11, 
respectively. 


MATCHING  OF  FIRST  TWO  MOMENTS 

We  now  derive  an  effective  number,  M^,  of  signal  bins,  in 
addition  to  an  effective  signal  power,  for  the  set  (S^)  for 
1  ^  m  5  M.  These  effective  parameters  would  then  be  used  in 
place  of  M  and  S  in  the  results  of  previous  sections,  and  would 
be  expected  to  have  wider  applicability  than  using  just  S^. 

For  nonlinearity  g(x)  ■  x  -  1  -  ln(x)  for  x  >  x^^,  define  the 
following  functions: 


a^{S,x^)  »  J  du  a  exp(-au)  [u  -  1  -  ln(u)]^  -  fj* 


(S.Xq) 


(69) 


“  tJ{S,x^)  -  #i(0,Xjj) 


(70) 
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Plots  of  the  three  functions  in  (70)  are  given  below,  nanely  the 
denominator,  the  numerator,  and  the  ratio  R(S,x^),  in  figures  60, 
61,  and  62,  respectively.  (If  desired,  the  integral  in  (69)  can 
be  evaluated  by  use  of  the  results  in  appendix  F.) 

The  mean  and  variance  of  the  output,  y^,  of  nonlinearity, 
g(x),  for  the  m-th  signal  input,  are  o  (S^^.x^) » 

respectively.  The  corresponding  first  two  cumulants  of  summer 
output  z  are 

M 

/ij  -  IZ  +  (N  -  M)  /i(0,x^)  - 

m*l 

My  * 

m*l 

-  C  *  **  ,  (72) 

respectively . 

On  the  other  hand,  if  all  the  signal  powers,  (S|||),  are 
replaced  by  an  effective  value,  for  an  effective  number,  M^, 
of  signal  bins,  we  obtain 

“  **e 

"z  “  **e  (‘^^<^'*0^  ■ 

Equating  (73)  to  (71),  and  (74)  to  (72),  it  follows  that  we  must 
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satisfy 

R{S,,X  )  .  - -  .  (75) 

IZ  (»'(S.«*o>  ■ 


Since  the  value  of  the  right-hand  side  of  (75)  is  )cnown,  once  M, 
{S|q)/  and  are  specified,  this  is  a  transcendental  equation  for 
effective  signal  power  That  is,  S  in  (70)  is  varied  until 

the  value  in  (75)  is  realized;  this  can  be  easily  done  by  use  of 
figure  62.  Then,  effective  number  follows  according  to 


It  should  be  noted  that  this  approach  accounts  for  the 
specific  nonlinearity  g(x),  the  breakpoint  value  x^,  and  all  of 
the  individual  signal  powers  for  1  6  m  1  M,  through  the 

quantities  in  (75),  which  in  turn  utilize  (63)  and  (69).  This 
approach  does  not  employ  an  ad  hoc  definition  of  effective 
values;  however,  it  does  require  a  numerical  solution  to  the 
transcendental  equation  (75)  for  S^,  or  the  use  of  figure  62. 

The  number  of  search  bins  M  is  not  altered  by  any  of  these 
replacements,  and  need  not  be  changed  in  utilizing  the  earlier 
plots.  Rather,  the  values  of  M  and/or  S  to  be  employed  in  using 
the  previously  plotted  receiver  operating  characteristics  are 
those  given  by  and/or  ^  here,  respectively. 
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SUMMARY 

The  M6LR  processor  has  been  analyzed  and  nusierically 
evaluated  in  terms  of  its  false  alarm  and  detection  probabilities 
Pf,Pd#  over  a  wide  range  of  parameter  values.  In  particular, 
plots  of  P^  versus  threshold  v  for  search  sizes  N  *  1024,  512, 
256,  128,  64  and  breakpoint  z^  «  1(1)11  have  been  presented. 

These  numerical  results  allow  for  accurate  determination  of 
required  threshold  settings  in  order  to  realize  specified  false 
alarm  probability  values. 

In  addition,  receiver  operating  characteristics,  P^  versus 
P£,  for  M  «  1024,  signal  set  sizes  M  •  8,  16,  32,  64,  128,  256, 
and  breakpoint  values  1«  3,  5,  7,  9,  11  have  been  plotted, 

with  the  signal  power  S  varied  over  a  range  sufficient  to  cover 
the  practical  range  of  performance  values.  All  of  these  results 
have  been  checked  by  simulation;  the  discrepancy  is  usually  less 
then  0.1  dB  in  terms  of  required  signal  power. 

The  energy  detector  has  been  analyzed  and  evaluated  exactly, 
and  can  serve  as  a  basis  of  comparison  with  the  M6LR  processor. 
Receiver  operating  characteristics  for  N  >  1024  and  M  >  8,  16, 

32,  64,  128,  256,  512,  1024  have  been  plotted.  The  corresponding 
plot  of  P£  versus  threshold  v  is  also  given. 

The  analysis  of  performance  has  been  accomplished  by  a 
combination  of  two  methods.  A  piecewise-linear  approximation 
works  best  for  large  breakpoint  values,  x^,  while  a  fourth-order 
moment  fitting  procedure  is  more  appropriate  for  small  x^  values. 
The  piecewise-linear  approach  can  actually  be  used  for  small 
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breakpoint  values,  such  as  ■■  3;  however,  the  tisw  of  execution 
is  sometimes  excessive,  and  there  is  the  danger  of  loss  of 
significance,  due  to  the  large  alternating  tenss  encountered  in 
the  series  expansions.  Accordingly,  the  fourth-order  moment  fit 
has  been  used  in  the  region  of  small  breakpoint  values  x^.  An 
alternative  is  investigated  in  appendix  6. 

One  possible  objection  to  use  of  these  results  for  the 
performance  of  the  M6LR  processor,  in  practice,  is  that  M  will 
often  not  be  kno«m,  meaning  that  the  best  breakpoint  x^  cannot  be 
selected.  This  raises  the  question  of  how  to  process,  in 
general,  when  M  is  unknown.  Also,  the  use  of  large  breakpoint 
values,  Xq,  (as  suggested  here  for  small  values  of  M)  causes  us 
to  reject  some  very  strong  signal  power  estimates,  which  is 
unsatisfactory  and  counter  intuitive.  In  fact,  the  question  of 
whether  we  should  employ  a  breakpoint  at  all  is  a  meaningful  one 
that  requires  further  investigation. 

It  is  worthwhile  to  derive  the  true  generalized  likelihood 
ratio  tests  for  the  situations  where  M  is  unknown,  as  well  as 
when  M  is  known  to  be  less  than  N,  but  the  signal  occupancy 
pattern  is  unknown.  Also,  if  partial  information  about  the 
signal  structure  is  available,  like  its  bandwidth,  how  should 
this  information  be  incorporated  in  the  generalized  likelihood 
ratio  test?  Comparisons  with  the  optimum  likelihood  ratio 
processors  or  near-optimum  processors  in  each  case  would  be  very 
worthwhile  and  should  be  pursued.  This  is  the  subject  of  a 
current  investigation  which  will  appear  in  a  NUWC  technical 
report  [ 3 ] . 
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APPENDIX  A.  MODIFIED  GENERALIZED  LIKELIHOOD  RATIO  PROCESSOR 

In  this  appendix,  we  will  derive  the  foxv  of  the  aodified 
generalized  likelihood  ratio  processor.  Since  the  actual  signal 
power  per  bin,  is  unkno%m,  we  hypothesize  paraaeter  value  S^ 
and  estimate  it  from  the  observation  set  (random  variables)  (x^). 
We  also  define  the  signal  power  parameter  a^^  ■  1/(1 
we  will  only  allow  signal  power  S^  1  0,  then  we  must  have  a^  1  1. 

According  to  (2),  the  joint  probability  density  function  of 
observation  (z^)  under  hypothesis  is,  for  hypothesized 
parameter  values 

N 

qj(Uj,...,Ujj)  -  TT  ®n  ^  ® 

n*l 

where  we  have  allowed  the  possibility  of  signal  being  present  in 
all  N  search  bins.  This  is  tantamount  to  assuming  N  «  N,  for  the 
time  being;  the  method  by  which  the  estimated  number  of  signal 
bins  is  reduced  below  the  value  N  will  be  indicated  below. 

For  observation  {x^},  the  density  in  (A-1)  takes  on  value 

N 

qj(Xj,...,Xj,)  -  TT  ®*P*"*n*n^*  • 

n*l 

This  quantity  is  to  be  maximized  by  choice  of  all  (Sq)*  It  can 
be  seen  that  the  best  choice  of  each  parameter  is  (random 
variable)  a^^  «  min { 1 , l/x^^ | ,  where  we  have  accounted  for  the 
constraint  a^  ^  1.  This  translates  into  tne  signal  power 
estimate  -  max|0,Xjj-it  -  However,  since  values  of  randan 
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variable  near  1  are  typical  (even  for  the  noise-alone  case  Bq 
where  B{Zj^)  is  1),  we  will  reject  signal  ponfer  estisates  which 
are  below  some  minimum  acceptable  level,  S^.  Thus,  our  siodified 
signal  power  estimate  is 


1  for  z_  -  1  2  S_ 
n  o 

0  otherwise 


Sq  i  0  . 


(A-3) 


This  is  equivalent  to  the  modification 


■  r: 


for  2  x_ 
n  o 

otherwise 


■  ; 


1  2  1 
o  o 


(A-4) 


Now,  let  L  be  the  set  of  n  values  in  (A-4)  for  which  z^  2  z^. 
Then,  upon  maximization,  the  n-th  term  in  (A-2)  becomes  l/{x^e) 
for  n  e  L,  and  ezp(-z^)  otherwise.  The  generalized  likelihood 
ratio  is  then  given  by 


6LR 


max  q.  _ 

—z -  -  TTlj-  exp(z  -1)1  -  TTlexptx  -l-ln(x  )])  - 

^o  nzL^^n  ”  *  ncL 

-  exp[j^[x^-l-ln(Xjj)])  -  ®*p(^  9(*„))  f  (A-5) 


where  nonlinearity  g(z)  is  defined  as 


g(x) 


1  -  ln(x) 
0 


for  X  2  Xj, 
otherwise 


(A-6) 


This  nonlinearity  has  a  breakpoint  at  x  ■  x^,  which  is  never  less 
than  1;  that  is,  from  (A-4),  x^  •  1  ♦  2  1,  since  minimum 

acceptable  signal  power  level  is  never  negative. 
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Coaparison  of  the  generalized  likelihood  ratio  in  (A>5)  with 
a  threshold  is  equivalent  to  the  test 


g(*n) 


(A-7) 


Since  breakpoint  of  nonlinearity  9(z)  rejects  any  input  values 

which  are  less  than  x^,  the  nusiber  of  nonzero  contributors  to 
n  o 

output  s  in  (A-7)  is  less  than  N.  In  fact»  if  is  significantly 
larger  than  1,  then  very  few  terois  will  contribute  to  the  output 
z  in  (A-7). 

The  selection  of  breakpoint  z^  in  this  M6LR  processor  is  its 
way  of  restricting  the  estimated  number  of  signal  bins  to  be  less 
than  M.  However,  unfortunately,  there  is  no  simple  way  of 
deciding  on  the  best  value  of  breakpoint  z^  for  each  hypothesized 
value  of  M.  We  must  conduct  a  complete  analysis  of  performance 
of  processor  (A-7)  in  terms  of  the  false  alarm  and  detection 
probabilities  in  order  that  the  best  z^  for  each  specified  N  can 
be  determined.  However,  this  still  does  not  address  the  problem 
of  what  to  do  when  M  is  unknown  in  practice.  This  is  a 
fundamental  limitation  of  the  M6LR  processor. 
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APPENDIX  B.  PIBCEWISB'LINBAR  FIT  TO  NONLINEARITY  9(z) 

Nonlinearity  9(z)  is  zero  for  z  <  z^,  and  otherwise  arbitrary 
for  z  2  z^,  where  z^  is  called  the  breakpoint.  Function  9(z)  is 
to  be  approzinated  by  a  piecewise-linear  fit  h(z),  which  is  also 
zero  for  z  <  z^,  but  which  is  linear  for  z  2  z^.  In  particular# 
let 

h(z)  -  c  +  b(z  -  z^)  for  z  2  z^  ,  (B-1) 

where  c  is  the  ordinate  of  h(z)  at  breakpoint  z^#  and  b  is  the 
slope  of  h(z)  for  z  >  z^.  A  variety  of  choices  for  b  and  c  are 
possible  and  will  be  developed  below. 

The  input  to  nonlinearity  9(z)  is  a  randosi  variable  x,  with 
probability  density  function  p^Cu).  The  (randosi)  error  between 
the  outputs  of  the  two  devices  is 

e  -  9(z)  -  h(z)  •  9(x)  -  c  -  b(z  -  x^)  for  z  2  z^  #  (B-2) 

and  e  >  0  for  z  <  z^.  Initially,  the  constants  b  and  c  in  fit 
h(z)  will  be  chosen  so  as  to  Binisdse  the  Bean  squared  error 

E  B  -  (g(z)  -  h(z)J^  *  J  P**®)  (9(“)  ”  *»(«))^  • 

a» 

-  J  du  Pjj(u)  (gCu)  -  c  -  b(u  - 
*o 

this  is  called  FIT-2,  because  two  paraBeters  are  optiBized. 

When  we  set  the  two  partial  derivatives  of  E  to  zero,  the 
best  values  of  b  and  c  are  found  to  be  given  by 
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j.  ^0  *0  •  *0  ^2  " 

b  - X —  ,  e  - - — 


lo  I-J  -  LJ 


Lft  L-  -  LT 
0  2  1 


(B-4) 


where 


-  J  du  Pjj{u)  (tt  -  , 


Rn  ■  J  P,(u)  g(u)  (u  -  . 


(B-5) 


In  the  special  case  where  probability  density  function  pj^(u) 
is  exponential  with  signal  power  paraaeter  a,  nasiely 


P«(u)  ■  a  exp(-au)  for  u  >  0  ,  a  ■ 


1 _ 

+  S  ' 


the  expression  for  simplifies  (for  1  0)  to 

-  nl  exp(-ax^)/a”  . 


(B-6) 


(B-7) 


Correspondingly,  the  optimum  parasieter  values  sis^lify  to 


B  -  a^  J  dt  exp(-at)  (at  -  1)  g(x^  +  t)  , 


S  -  a  J  dt  exp(-at)  (2  -  at)  g(xQ  +  t)  . 


(B-8) 


These  optimum  values,  b  and  S,  obviously  depend  on  both  the 
breakpoint  x^  as  well  as  the  probability  density  signal  power 
parameter  a.  For  signal  power  S  «  0,  then  a  ■  1,  and  B,  S  are 
denoted  by  B^,  S^,  which  depend  only  on  x^. 
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There  are  cases  where  it  is  desirable  to  fix  ordinate  c,  a 
priori,  and  to  nininize  mean  square  error  B  solely  by  choice  of 
slope  b;  this  is  called  FIT-1.  For  ezanple,  ordinate  C  above  can 
turn  out  to  be  negative,  even  though  original  nonlinearity  g(z) 
is  always  positive  for  z  2  z^.  If  this  is  undesirable,  c  could 
be  set  to  zero,  for  example,  and  the  best  slope  b  selected.  When 
we  set  the  derivative  of  mean  square  error  E  (with  respect  to  b) 
equal  to  zero,  the  optimum  slope  (for  fixed  c  ■  c)  is  given  by 


E  -  J  du  Pjj(u)  (g(u)  -  c]  (u  -  x^Y  J  du  p^(u)  (u  - 
z_  '  z_ 


*o>' 


(B-9) 


In  the  special  case  of  an  exponential  probability  density 
function  p^(u) ,  as  in  (8-6)  above,  there  follows  (for  x^  2  0) 


B  ■  y  a^  J  dt  exp(-at)  t  (g(XQ  +  t)  -  c]  .  (B-10) 

0 

This  optimum  value,  E,  depends  on  x^,  a,  and  c.  For  sigxml 
power  S  "  0,  then  a  *  1,  and  we  denote  c  by  c^,  and  E  by  E^ 
which  depends  on  x^  and  c^. 

Finally,  the  parameters  b  and  c  in  fit  h(x)  can  be  chosen 
arbitrarily,  with  no  specific  optimization  in  mind;  this  case, 
where  zero  parameters  are  optimized,  is  called  FIT-0.  For 
example,  we  might  take  c  >  9(Xq)#  the  initial  ordinate,  and  take 
b  «  g'(x^),  the  initial  slope;  this  corresponds  to  the  tangent 
line  to  g(x)  at  x^.  (This  particular  choice  takes  no  account  of 
the  signal  power  S. )  Generally,  b  and  c  can  depend  on  parameter 
a  as  well  as  x^.  When  S  ■  0,  then  a  -  1,  and  we  denote  b  and  c 
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by  and  c^,  respectively,  which  depend  only  on  z^.  Sooe  of  the 
coefficients,  c  and  b,  are  tabulated  below  for  the  various  fit 
procedures . 

The  mean  square  error  B  for  FIT-2  is  always  smaller  than 
those  for  the  other  two  fits.  However,  there  is  the  possibility 
of  a  negative  ordinate  &  or  S^,  which  may  or  may  not  be  damaging. 
On  the  other  hand,  FIT-1  may  be  poor  for  large  signal  powers  S, 
where  large  values  of  random  variable  z  are  frequently 
encountered.  Finally,  FIT-0  can  be  very  unsatisfactory,  (as  for 
example  by  using  the  tangent,)  since  the  output  of  h(z)  might 
always  underestimate  the  true  nonlinearity  g(z);  this  could  be 
especially  damaging  for  large  S. 

One  way  around  this  shortcoming  of  the  approzimation  h(z) 
used  in  (B-1)  is  to  use  several  breakpoints  in  approzimating  g(z) 
with  piecewise-linear  fits.  One  possible  problem  with  this 
latter  approach  is  the  tractability  of  the  Fourier  transform  of 
input  probability  density  function  p^Cu)  over  an  arbitrary  finite 
segment.  For  an  exponential  probability  density  function,  this 
is  not  a  problem.  However,  a  more  significant  problem  will 
arise  in  determining  the  dominant  asymptotic  decay  of  the 
characteristic  function  of  the  output  of  nonlinearity  h(z),  so 
that  it  can  be  subtracted  out  for  rapid  convergence  of  the  FFT 
required  for  numerical  evaluation  of  the  exceedance  distribution 
function.  This  procedure  has  already  been  carried  out  for  a 
nonlinearity  with  a  single  breakpoint.  Although  tedious,  this  is 
a  workable  approach  for  accurately  approzimating  the  detection 
characteristics  of  any  nonlinear  memoryless  device. 
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OPTIMUM  COEFFICIENT  b  FOR  FIT-1,  g(x)  -  z 


s 

b 

C 

0 

1 

.50000  00000 

1 

1 

.63463  61709 

0 

2 

1 

.70479  93293 

0 

3 

1 

.74929  19895 

0 

4 

1 

.78053  21002 

0 

5 

1 

.80390  29111 

0 

6 

1 

.82216  51352 

0 

7 

1 

.83689  78455 

0 

8 

1 

.84907  68432 

0 

9 

1 

.85934  10855 

0 

10 

1 

.86812  82760 

0 

16 

1 

.90265  28405 

0 

0 

2 

.83409  07182 

0 

1 

2 

.82671  32049 

0 

2 

2 

.84136  46500 

0 

3 

2 

.85567  46879 

0 

4 

2 

.86781  56014 

0 

5 

2 

.87797  07329 

0 

6 

2 

.88652  82066 

0 

7 

2 

.89382  42960 

0 

8 

2 

.90011  88888 

0 

9 

2 

.90560  86911 

0 

10 

2 

.91044  31283 

0 

16 

2 

.93065  54450 

0 

0 

2 

.68066  43084 

9(X, 

1 

2 

.75000  00000 

9(x; 

2 

2 

.79022  25134 

9(x[ 

3 

2 

.81731  80855 

9(Xf 

4 

2 

.83713  03195 

9(x 

5 

2 

.85239  96646 

9(Xf 

6 

2 

.86461  01481 

9(x 

7 

2 

.87464  59948 

9(xJ 

8 

2 

.88307  15099 

9(xJ 

9 

2 

.89026  60501 

9(xJ 

10 

2 

.89649  52728 

9(xJ 

16 

2 

.92163  03621 

9(xJ 

-  1  -  ln(x) 


0 


l-ln(2) 
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OPTIMUM  COEFFICIENT  b  FOR  FIT-1,  g(x)  -  x  - 


X 


o 


b 


c 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 


50000 

00000 

68066 

43084 

76208 

37403 

80951 

84749 

84084 

43526 

86316 

90731 

87992 

32882 

89297 

87374 

90344 

78223 

9(*o 

91203 

50027 

91920 

88943 

1  -  ln(x) 
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OPTIMUM  COEFFICIENTS  FOR  FIT-2,  g(x)  -  x  -  1  -  ln(x) 


dB) 

*o 

—• 

11 

-4 

11 

-3 

11 

-2 

11 

-1 

11 

0 

11 

1 

11 

2 

11 

3 

11 

4 

11 

5 

11 

6 

11 

7 

11 

8 

11 

9 

11 

10 

11 

9 

-4 

9 

-3 

9 

-2 

9 

-1 

9 

0 

9 

1 

9 

2 

9 

3 

9 

4 

9 

5 

9 

6 

9 

7 

9 

8 

9 

9 

9 

10 

9 

E 

.92225  95675 
.92620  52679 
.92714  32728 
.92828  10207 
.92964  98384 
.93128  09423 
.93320  30816 
.93543  95228 
.93800  45899 
.94090  01605 
.94411  26969 
.94761  14546 
.95134  84273 
.95526  03258 
.95927  24991 
.96330  43057 

.90775  76002 
.91315  09831 
.91441  81751 
.91594  78247 
.91777  76590 
.91994  35729 
.92247  62583 
.92539  71647 
.92871  41971 
.93241  77882 
.93647  81265 
.94084  42937 
.94544  58154 
.95019  66890 
.95500  14422 
.95976  23429 


e 

7.5960  03381 
7.5912  90368 
7.5899  30861 
7.5881  49569 
7.5858  06201 
7.5827  13409 
7.5786  23271 
7.5732  12010 
7.5660  64105 
7.5566  57775 
7.5443  54565 
7.5283  96232 
7.5079  11992 
7.4819  38230 
7.4494  51068 
7.4094  10003 

5.7941  55867 
5.7877  17092 
5.7858  80522 
5.7834  85744 
5.7803  53285 
5.7762  46689 
5.7708  57817 
5.7637  91177 
5.7545  49005 
5.7425  19630 
5.7269  72270 
5.7070  61457 
5.6818  43574 
5.6503  06459 
5.6114  10871 
5.5641  40375 
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OPTIMUM  COEFFICIENTS  FOR  FIT-2,  g(x)  -X  -  1  -  ln(x) 


S((iB) 

*o 

— «• 

7 

.88648 

-7 

7 

.89059 

-6 

7 

.89160 

-5 

7 

.89283 

-4 

7 

.89433 

-3 

7 

.89614 

-2 

7 

.89832 

-1 

7 

.90090 

0 

7 

.90393 

1 

7 

.90744 

2 

7 

.91144 

3 

7 

.91592 

4 

7 

.92086 

5 

7 

.92619 

6 

7 

.93184 

7 

7 

.93769 

8 

7 

.94363 

9 

7 

.94954 

10 

7 

.95529 

•.go 

5 

.85211 

-4 

5 

.86465 

-3 

5 

.86748 

-2 

5 

.87085 

-1 

5 

.87481 

0 

5 

.87940 

1 

5 

.88464 

2 

5 

.89052 

3 

5 

.89700 

4 

5 

.90400 

5 

5 

.91142 

6 

5 

.91910 

7 

5 

.92690 

8 

5 

.93465 

9 

5 

.94219 

10 

5 

.94938 

C 


76725 

4.0409 

61082 

83768 

4.0364 

46787 

15654 

4.0352 

17624 

16788 

4.0336 

38712 

08847 

4.0316 

04404 

47388 

4.0289 

75636 

03044 

4.0255 

69793 

32451 

4.0211 

48392 

38303 

4.0154 

02859 

20032 

4.0079 

39227 

19462 

3.9982 

63354 

68934 

3.9857 

69096 

51677 

3.9697 

32518 

84553 

3.9493 

15412 

30742 

3.9235 

80802 

44685 

3.8915 

21627 

44925 

3.8521 

01661 

04350 

3.8043 

05370 

43858 

3.7471 

91510 

08814 

2.3680 

29030 

39431 

2.3530 

81233 

93493 

2.3489 

72161 

95380 

2.3436 

96494 

86148 

2.3369 

20078 

72956 

2.3282 

21660 

55460 

2.3170 

78963 

54537 

2.3028 

58079 

56821 

2.2848 

09313 

89311 

2.2620 

72753 

35093 

2.2336 

96167 

93940 

2.1986 

66239 

81996 

2.1559 

51865 

55879 

2.1045 

55788 

41367 

2.0435 

68946 

46865 

1.9722 

21122 

114 


TR  10539 


OPTIMUM  COEFFICIENTS  FOR  FIT-2,  g(x)  ■  X  -  1  -  ln(x) 


s(dB)  Xj,  E  e 


3 

.78625 

12208 

.85305 

27503 

-4 

3 

.80999 

00501 

.82481 

06427 

-3 

3 

.81516 

25120 

.81731 

52922 

-2 

3 

.82122 

36047 

.80782 

83112 

-1 

3 

.82822 

71199 

.79584 

26744 

0 

3 

.83619 

25020 

.78074 

60381 

1 

3 

.84509 

41622 

.76181 

49196 

2 

3 

.85485 

37960 

.73821 

72304 

3 

3 

.86533 

79676 

.70902 

61829 

4 

3 

.87636 

25481 

.67324 

74273 

5 

3 

.88770 

42542 

.62985 

91859 

6 

3 

.89911 

79023 

.57786 

25469 

7 

3 

.91035 

66214 

.51633 

66733 

8 

3 

.92119 

16319 

.44449 

22468 

9 

3 

.93142 

85042 

.36171 

64722 

10 

3 

.94091 

78565 

.26760 

44129 

11 

3 

.94955 

98168 

.16197 

38256 

12 

3 

.95730 

28205 

.04486 

33282 

13 

3 

.96413 

81279 

-.08348 

39910 

14 

3 

.97009 

17657 

-.22264 

79238 

15 

3 

.97521 

54947 

-.37206 

52497 

1 

.59634 

73623 

-.19269 

47246 

-4 

1 

.66461 

29941 

-.27361 

79847 

-3 

1 

.67837 

13161 

-.29355 

18221 

-2 

1 

.69403 

34922 

-.31806 

14225 

-1 

1 

.71154 

20343 

-.34801 

67111 

0 

1 

.73072 

76581 

-.38436 

59487 

1 

1 

.75130 

37657 

-.42810 

59470 

2 

1 

.77287 

66361 

-.48023 

95488 

3 

1 

.79497 

21515 

-.54172 

30324 

4 

1 

.81707 

57330 

-.61340 

88469 

5 

1 

.83867 

85685 

-.69599 

10057 

6 

1 

.85932 

16212 

-.78996 

07593 

7 

1 

.87862 

98970 

-.89557 

85961 

8 

1 

.89633 

24642 

-1.0128 

65645 

9 

1 

.91226 

74477 

-1.1416 

14117 

10 

1 

.92637 

43448 

-1.2814 

13515 
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APPENDIX  C.  PARTIAL  FRACTION  EXPANSION  OF  PRODUCT 


For  j  and  n  integer,  and  coo^lex  paraaeters  a  fi,  define  the 
quantities 


Vjn(a,P) 


(a-p) 


(fi-a) 


<J>n 

tin 

n! 

<J'n 

3+n 

n! 

for  j  2  1  ,  n  2  0  ; 


for  j  2  1  ,  n  2  0 


(C-1) 


(C-2) 


Then,  »jn(P»a)  -  #ij„(a,p). 

These  quantities  can  be  easily  evaluated  by  recursions: 


(o^)^  for  j  i  1  , 

p~a  n 

.jil; 

(C-3) 

Vjg(a,|l)  - 

(Af  to"-  3  »  1  ' 

«  for  n  i  1 

a-p  n 

,  j  i  1  . 

(C-4) 

Now,  consider  the  partial  fraction  expansion  of  the  product 
of  the  following  two  functions,  where  complex  constants  a  S; 
for  integers  m2  1,  k2  1,  define  the  function  of  complex 
variable  z. 


f(z)  ■  —  •  (C-5) 

(l+za)“  (1+zS)* 

Develop  f(z)  according  to  partial  fraction  expansion 
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m-1  n  k-1 

«(*)  -  El - ^  ’ 

n»6  (l+za)"  “  n«d  ( 


(C-6) 


Multiply  both  sides  of  (C>6)  by  (l**‘za)*  to  obtain 


^  St"!  _ 

- - — n  -  LZ  «n  ♦  (l+so)“ 

(l+z(l)*  n«0  ”  n*6 


k-l 

C - h=E  • 

(l+xp) 


Take  the  p-th  derivative  of  (C-7)  with  respect  to  z,  and  then  let 
z  -  -1/a;  the  right-hand  side  becomes  simply  Up  pi  while  the 
left-hand  side  is 


(-»)“  ,K, 


(C-8) 


Equating  the  two  sides,  we  recognize  that  coefficient  Up  in  (C-6) 
is  equal  to  A>}(p(a«6)  defined  in  (C-1).  A  similar  treatment  for 
coefficients  (v^)  in  (C-6)  yields  the  following  expansion:  for 
m  2  1,  k  1  1,  a  6, 


(1+za)®  (l+zP)^  ^  (1+zo)®"”  n=5  (1+zP)*^"" 


v^(o,p) 

-)iZn  • 
1+zP)*  “ 


Observe  that  expansion  coefficients  6) )  depend 

on  m,  and  that  expansion  coefficients  {Vjg2^(a,6))  do  not  depend  on 
k.  It  should  also  be  noted  that  numerical  difficulties  can  be 
expected  with  series  expansion  (C-9)  when  m  and  k  are  large,  and 
o  is  close  to  6,  due  to  the  alternating  character  of  the  terms. 
Although  the  terms  themselves  can  be  calculated  very  accurately 
by  means  of  (C-1)  and  (C-2),  the  sums  in  (C-9)  will  lose 
significance  due  to  cancellation  effects. 
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APPENDIX  D.  TWO  ALTERNATIVE  EXACT  ANALYSIS  NETBOOS 
HIGH-ORDER  MOMENT  APPROACH 

We  can  use  the  methods  of  (4;  page  88]  to  determine  the  exact 
exceedance  distribution  function  of  sussMr  output  s  in  (4), 
provided  that  we  can  accurately  evaluate  the  higher-order  moswnts 
of  the  summer  inputs  These  mcxsents  can  be  obtained  as 

follows;  for  x^  1  1,  and  g(x)  •  x  -  1  -  ln(x)  for  x  >  x^,  zero 
otherwise,  the  k-th  moment  of  is  given  exactly  by 

_  m 

^n  "  -  Jdu  Pjj(u)  g(u)*‘  •  J  du  Pjj(u)  (u  -  1  -  ln(u))*^  (D-1) 

*o 

for  k  1  1.  For  the  exponential  density  (2)  of  input  random 
variables  this  k-th  moment  becoses 

yj^  «  J  du  a  exp(-au)  (u  -  1  -  ln(u))^  for  k  k  1  ,  (D-2) 

*o 

which  depends  on  k,  x^,  and  a.  Since  the  integrand  is  positive 
for  all  u,  regardless  of  k  and  a,  these  integrals  can  be 
conducted  very  accurately  to  high  order  k,  and  used  in  the 
recursive  procedures  of  [4]  to  get  the  cusnilants  of  {y^}  and  s. 

The  major  stumbling  block  to  this  approach  is  the  round-off 
error  encountered  in  transforming  moments  to  cumulants.  Double 
precision  arithmetic  is  a  necessity  for  useful  numerical  results. 
For  k  «  1,  a  closed  form  result  is  available;  see  (63).  For 
k  «  2,  (D-2)  could  be  evaluated  from  the  results  in  appendix  F. 


119 


TR  10539 


CHARACTERISTIC  FUNCTION  APPROACH 

An  alternative  approach  is  to  directly  nuaerically  evaluate 
the  characteristic  function  of  lyj^l#  which  is  given  exactly  by 

•  exp(i^yjj)  -  exp(i5g(Xjj))  ■  J  du  p^^Cu)  exp(iCg(u))  - 
*o 

■  J  du  Pjj(u)  +  J  du  Pjj(u)  exp(iC[u  -  1  -  ln(u)])  .  (D-3) 

*o 

When  probability  density  function  p^Cu)  is  exponential,  this 

reduces,  for  x^  >  0,  to 
o 

«D 

fy(C)  »  1  -  expC-ax^)  +  a  J  du  exp(-au  +  it[u  -  1  -  ln(u) ]) . (D-4) 

*o 

The  integral  depends  on  K,  a,  and  x^.  It  snist  be  nuaerically 
evaluated  very  accurately  for  nuaerous  values  of 

The  aajor  drawback  with  this  approach  is  the  lack  of  decay  of 
characteristic  function  fy(()  with  K*  Nasiely,  fy(C)  behaves  as 
the  constant,  1  -  exp(-ax^),  as  ^  ±«».  In  order  to  nuaerically 

Fourier  transform  f-(^),  which  is  given  in  teras  of  powers  of 
f  must  determine  the  dominant  asymptotic  behaviors  of 
f  (C)  auid  f-(C)-  In  fact,  it  will  be  necessary  to  determine 
several  of  the  leading  terms  in  the  asymptotic  expansion  of 
C ) f  least  to  the  point  where  the  remainder  decays  fast 
enough  to  ensure  that  the  termination  of  the  Fourier  integral,  on 
the  remainder  of  f-(^),  contributes  insignificantly  to  the  trun- 
cation  error.  The  asymptotic  behavior  of  fy(C)  is  now  considered. 
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ASYMPTOTIC  BEHAVIOR  OF  CHARACTERISTIC  FUNCTION 
FOR  IMPULSIVE  PROBABILITY  DENSITY  FUNCTION 

The  generic  problem  is  es  follows x  a  given  probability 
density  function  Py(u)  has  an  isipulse  of  area  A^  at  the  origin, 
and  a  distributed  cosiponent  which  is  nonzero  only  for  u  1  y^, 
where  y^  -  g(*Q).  Furthenoore,  Py(u)  is  smooth  for  u  >  y^  and 
decays  to  zero  for  large  u*  The  corresponding  characteristic 
function  is 

«• 

tyil)  -  Aj  ♦  J  du  ezp(Uu)  py(u)  .  (D-5) 

^o 

We  integrate  by  parts,  letting  (for  u  >  y^) 

U  *  exp(oni)  Py(u)  ■  h(u)  ,  dV  •  du  ezp(-(a>il)tt)  ,  (D-6) 

where  positive  real  constant  o  is  yet  to  be  chosen.  Then, 

du  -  du  h'(u)  ,  V  -  -  ,  (D-7) 

leading  to 

exp(-(a-i5)y_)  ,  7 

fy(^)  •  - 2-  h(y^)  +  Jdu  exp(-(o-U)u)  h'(u). 

^o  (D-8) 

Repeated  integrations  by  parts  yields  the  asymptotic  expansion  of 
characteristic  function  fy{K)  as 

f„(^)  ~  A,  +  El - exp  (-ay  )  h»  ^yj,)  as  ^  ±»  . 

y  ^  n-1  (a-U)"  ”  (D-9) 

Strictly,  the  functions  and  derivatives  are  evaluated  at  u  •  y^^. 
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Becaus*  b(u)  •  azp(cQi)  Py(tt),  the  darivativaa  raguirad  above 
can  ba  avaluatad  as  follows;  for  u  >  y^, 

h'(tt)  •  azp(otu)  (o  Py(tt)  ♦  Py(u)J  ,  (D-10) 

b*(u)  -  exp(au)  [a^  Py(tt)  ♦  2o  pj(u)  ♦  pj(tt)]  .  (D-11) 

Then,  the  asyaptotic  expansion  for  characteristic  function  fy(C) 
becoses  explicitly,  to  third  order,  f(()  f3(()  as  C  idiere 


Aj  ♦  «p(«y,)  ^ 


il/ar  (1-U/o) 


/  (D-12) 


Pj  -  5  Py(yo)  *  ^2  "  4(®  Py<yo>  Pi^^oO  ' 


^3  •  Py^^o*  **“  Py<yo>  **“  Py<yo>)  * 


One  possible  choice  for  a  is  to  aake  the  coefficient,  ^2^  of 
—2 

(l-i5/o)  equal  to  zero;  that  is,  take  a  ■  ”Py(yo)/Py(yo) »  or 
approximately  so.  Another  is  to  let  exp  (-cm)  duplicate  the 
asymptotic  behavior  of  Py(u);  that  is,  have  h(u)  approach  a 
constant  as  u  'f*.  The  choice  of  a  does  not  seem  critical. 

The  individual  terms  in  f3(^)  can  be  easily  Fourier 
transformed  analytically.  The  leading  term  is  an  impulse  at 
u  •  0,  as  before.  The  next  tens  is  a  decaying  exponential 
starting  at  u  ■  y^.  The  (l-ill/a)”^  term  leads  to  a  u  exp(-u) 
type  term,  also  delayed  to  u  «  y^.  The  last  term  behaves  as 
u  exp(-u),  but  starting  at  u  ■  y^. 
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We  are  also  interested  in  the  asyaptotic  decay  ot  a  po«ier  of 
Letting  z  ■  (l-iC/a)“^,  we  have  the  k-th  power 

Ic 

«y(5)  -  -  K  ♦  «xp(Uyj,){Pi  Z  ♦  Pj  z^  ♦  P3  z^))  - 

-  aJ  +  k  aJ‘^  exp(iCyj,|(Pj  z  ♦  Pj  z^  +  P3  z^)  ♦  (J)  aJ‘^  x 

X  exp(iUyjj)(pi  2  Pj  Pj  *^)  ♦  (5)  exp(iC3y^)  pj  z^  , 

(D-15) 


as  z  *»  0,  to  third  order.  Substituting  for  z. 


tyil)  -  f3,^(C)  ■  aJ  +  k  aJ”^  exp(Uy^) 


T=IUa  * 


(l-i^/o) 


5 


(1-U/o) 


+  (3I  exp(iUy-)  - ^  as  5  -►  ±-  .  (D-16) 

This  is  the  desired  result  for  the  asysg>totic  behavior  of  the 
k-th  power  of  characteristic  function  fy(C)> 

The  corresponding  cumulative  distribution  function  is 

aJ  Cq(v)  +  k  aJ“^(pj  Cj(Vj)  +  P2  CjCVj)  +  P3  CjCVj))  + 

♦  (5)  aJ'"  (p?  C2(Vj)  +  2  Pi  Pj  C3(Vj))  +  ($)  aJ'*  pJ  CjCVj)  , 

(D-17) 

where 

»  a(v  -  ny_)  .  (D-18) 

n  o 
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EVALUATION  OF  DERIVATIVES  OF  Py(U) 

In  order  to  use  the  results  in  (D>12)  and  (D-14),  we  need  to 

know  the  derivatives  of  probability  density  function  Py(u)  at 

u  ■  yo”*"'  Yq  "  9^*0^*  Randon  variable  y^^  is  obtained  by 

nonlinear  transformation  where  nonlinearity  9(x)  is 

presumed  to  be  strictly  monotonically  increasing  for  z  >  and 

to  be  zero  for  z  <  z^.  We  now  relate  the  desired  derivatives  of 

o 

Py(u)  to  those  of  Pj^(u)  and  the  nonlinear  transformation  9(z). 
This  will  enable  us  to  determine  the  constants  F3 

required  in  (D-14). 

For  V  >  the  exceedance  distribution  function  of  random 
variable  y^  is  given  by 

E  (V)  -  Pr(y^  >  V)  -  Pr(g(z^)  >  v)  •  >  9(v))  -  Ejj(9(v))  . 

(D-19) 

Here,  g(y)  is  the  inverse  function  to  g(x),  at  least  for  y  >  y^; 
that  is, 

if  y  -  g(z)  for  z  >  z^  ,  then  z  -  9(y)  for  y  >  y^  -  gCz^)  . 

(D-20) 

A  useful  relation  follows  from  the  fact  that  g(§(y))  ■  y*  Taking 
a  derivative  with  respect  to  y,  we  have  g'(9(y))  9'(y)  ■  1  ,  or 

-iMlIyTT  • 

The  probability  density  function  of  random  variable  y^^  is 
given  by 
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Py(U) 


^  By(u)  -  -  ^  B^(5(u))  -  Pjj(9(u))  9'(U)  - 


Pje(9(u)) 

"  g*(^(u)) 


for  u  >  , 


(D-22) 


where  we  used  (D-21).  If  we  now  let  u  «  9(t},  this  latter 
relation  takes  on  the  desired  very  useful  fora: 


P,(t) 

Py(g{t))  ■  •  (D-23) 

Substitution  of  t  ■  yields  PyCy^)  “  ^*o^ '  which  is 

one  of  the  results  required  for  (D-14).  Notice  that  inverse 
function  §(y)  is  not  required  or  used  in  (D>23)  or  the  sequel. 

We  now  take  two  successive  derivatives  of  (D-23)  with  respect 
to  t,  to  get 


P^(g(t)) 


g*(t)  P4(t)  -  p^(t)  g"(t) 


(D»24) 


Pj(g(t)) 


3  g'  g"  p; 


3  g"‘  -  g'  g' 


i 


,.5 


(D-25) 


where  all  the  quantities  in  the  last  equation  are  to  be  evaluated 
at  argument  t.  When  we  substitute  t  ■■  we  obtain  expressions 
for  Py(yQ)  and  Pj(yo^  terms  of  quantities  easily  available 
from  the  given  functions  Pjj(n)  and  g(x).  There  is  no  need  to 
evaluate  inverse  function  9(y)  or  its  derivatives  at  all. 

For  example,  we  have,  for  x  >  x^^  and  for  u  >  0, 
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g(x) 

P^(u)  - 


■  X  -  1  -  ln(x),  g'(x)  “  ^  ”  X  '  S"!*)  ■  "5  »  (D-26) 


exp(-au)  ,  p^(u)  ■  -a^  exp(-au)  ,  P^(u)  -  a^  exp(-au)  . 
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APPENDIX  B.  TABLES  OP  FIRST  FOUR  CUMULANTS  OF 


FIRST  FOUR  CUMULANTS  OF 

y„  roR 

M 

o 

1 

S(dB) 

Xy(2,Xo,S) 

Xy(3yXo'-* 

Xy(<» 

.14849 

55068 

.17363 

54789 

.36390 

41244 

1.1105 

77108 

-3 

.37232 

63923 

.67739 

77766 

2.1866 

94609 

10.222 

37915 

-2 

.44083 

61060 

.87450 

33444 

3.0742 

91858 

15.641 

93753 

-1 

.53192 

32228 

1.1641 

73519 

4.5105 

69661 

25.281 

35153 

0 

.65328 

77247 

1.5963 

40984 

6.9002 

68199 

43.135 

87531 

1 

.81510 

93535 

2.2499 

35901 

10.982 

78191 

77.532 

33382 

2 

1.0307 

05214 

3.2500 

85125 

18.129 

05704 

146.29 

59138 

3 

1.3173 

12381 

4.7954 

26841 

30.913 

39282 

288.51 

60247 

4 

1.6970 

24416 

7.2014 

17458 

54.221 

46053 

591.76 

61973 

5 

2.1979 

28595 

10.968 

93324 

97.408 

56056 

1255.9 

48815 

6 

2.8555 

11925 

16.892 

91290 

178.52 

20929 

2744.8 

30537 

7 

3.7144 

26320 

26.234 

42995 

332.59 

30788 

6149.3 

66839 

8 

4.8307 

25577 

40.993 

28950 

627.96 

80417 

14066. 

27845 

9 

6.2747 

10254 

64.339 

96462 

1198.5 

59409 

32737. 

77500 

10 

8.1345 

44258 

101.30 

01036 

2307.6 

61144 

77293. 

54256 

11 

10.520 

84245 

159.83 

93733 

4474.4 

34247 

18465  ! 

5.1198 

12 

13.572 

47552 

252.58 

28253 

8724.9 

42288 

44542  8.8106 

13 

17.463 

89993 

399.53 

99358 

17090. 

91588 

10829  1 

B1.677 

FIRST  FOUR  CUMULANTS  OF 

y„  tor 

Xq  -  2 

S(dB) 

Xy(l/; 

Xy(2,XQ,S) 

Xy(3/*o'-^ 

Xy(<y: 

.12796 

27858 

.17574 

47137 

.37421 

11513 

1.1387 

25648 

-3 

.34822 

40437 

.69043 

82072 

2.2305 

06064 

10.384 

05135 

-2 

.41655 

24216 

.89094 

29992 

3.1298 

99242 

15.860 

39913 

-1 

.50764 

21649 

1.1850 

11999 

4.5823 

63547 

25.584 

30823 

0 

.62925 

97917 

1.6227 

77524 

6.9941 

24512 

43.565 

01106 

1 

.79164 

09435 

2.2833 

42973 

11.106 

35405 

78.149 

58412 

2 

1.0081 

35006 

3.2919 

45145 

18.291 

97917 

147.19 

16125 

3 

1.2959 

77672 

4.8472 

32578 

31.127 

31394 

289.81 

88541 

4 

1.6772 

24228 

7.2645 

48258 

54.499 

84453 

593.65 

45829 

5 

2.1798 

91997 

11.044 

52901 

97.766 

23265 

1258.6 

63214 

6 

2.8393 

78415 

16.981 

77103 

178.97 

45759 

2748.6 

86124 

7 

3.7002 

41138 

26.336 

94483 

333.15 

58187 

6154.7 

66460 

8 

4.8184 

48437 

41.109 

44269 

628.65 

56484 

14073. 

72632 

9 

6.2642 

33059 

64.469 

36116 

1199.3 

85027 

32747. 

89234 

10 

8.1257 

12529 

101.44 

20408 

2308.6 

36021 

77307. 

08684 

11 

10.513 

47608 

159.99 

29258 

4475.5 

67530 

18467 

3.0099 

12 

13.566 

38612 

252.74 

69308 

8726.2 

41019 

44545 

2.1601 

13 

17.458 

90366 

399.71 

34719 

17092. 

38514 

10830 

11.838 
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FIRST  FOUR  CUMULANTS  OF  FOR  *  3 


S(dB) 

Xy(2,Xo,S) 

Xy(3»*o'S* 

Xy(<r 

.08161  61389 

.15797  37056 

.38690  96395 

1.2208 

51895 

-3 

.27666  02014 

.69156  35364 

2.3702  27986 

11.036 

89414 

-2 

.34121  53043 

.90194  63156 

3.3221  83512 

16.774 

78383 

-1 

.42874  07666 

1.2104  17015 

4.8501  33315 

26.896 

18335 

0 

.54736  75223 

1.6685  94431 

7.3691  87760 

45.481 

75146 

1 

.70773  91341 

2.3571  33468 

11.631  37379 

80.982 

51465 

2 

.92358  84667 

3.4026  30423 

19.021  84015 

151.39 

81523 

3 

1.2124  28857 

5.0045  81998 

32.129  31540 

296.05 

40638 

4 

1.5964  18859 

7.4783  77694 

55.852  29933 

602.83 

07406 

5 

2.1034  74853 

11.323  74546 

99.555  41783 

1272.0 

15276 

6 

2.7687  09493 

17.333  46922 

181.29  04235 

2767.8 

40967 

7 

3.6362  73368 

26.765  75330 

336.08  70650 

6181.8 

14926 

8 

4.7616  84723 

41.617  23609 

632.28  52476 

14111. 

29939 

9 

6.2147  55389 

65.055  34614 

1203.7  86484 

32799. 

24592 

10 

8.0832  59751 

102.10  31337 

2313.8  71014 

77376. 

20297 

11 

10.477  54291 

160.72  42975 

4481.6  84888 

18476 

4.7253 

12 

13.536  32086 

253.54  25985 

8733.2  76884 

44557 

2.3351 

13 

17.433  99077 

400.56  68463 

17100.  36401 

10831 

67.581 

FIRST  : 

POUR  CUMULANTS  OF 

FOR  •  4 

■'ll  o 

S(dB) 

Xy(l»Xo'^^ 

Xy(2/*Q»S) 

Xy(3»XQ,S) 

Xy(4/ 

.04409 

23362 

.11669  06448 

.35266  30056 

1.2490 

58964 

-3 

.19680 

45752 

.62991  50503 

2.4586  33259 

11.977 

52732 

-2 

.25275 

52920 

.84403  91675 

3.4874  09776 

18.216 

44483 

-1 

.33098 

77685 

1.1623  08487 

5.1381  14630 

29.139 

34394 

0 

.44016 

13655 

1.6398  93287 

7.8490  33261 

48.998 

37207 

1 

.59167 

89225 

2.3618  68137 

12.402  31198 

86.492 

11951 

2 

.80018 

53326 

3.4596  25997 

20.218  49273 

159.95 

79889 

3 

1.0841 

03120 

5.1369  50753 

33.922  84360 

309.15 

65917 

4 

1.4663 

01430 

7.7119  13587 

58.445  51834 

622.50 

32500 

5 

1.9750 

40621 

11.684  38786 

103.17  17742 

1300.9 

23790 

6 

2.6453 

34913 

17.844  38888 

186.15  89727 

2809.4 

09138 

7 

3.5208 

35504 

27.444  78839 

342.42  69507 

6240.3 

71327 

8 

4.6562 

93538 

42.475  26777 

640.29  29312 

14192. 

26152 

9 

6.1206 

54795 

66.095  85133 

1213.6  27511 

32909. 

34237 

10 

8.0008 

77841 

103.32  27879 

2325.6  76461 

77523. 

73434 

11 

10.406 

63981 

162.11  42758 

4495.5  51594 

18495 

9.8717 

12 

13.476 

17733 

255.09  01757 

8749.2  70217 

44582 

7.5304 

13 

17.383 

59317 

402.25  69708 

17118.  52208 

10834 

97.984 
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FIRST  FOUR  CUMULAMTS  OF 

y^FOR 

*o  •  5 

S(dB) 

Xy( ^ 

Xy(2,Ji 

Xy(3#a 

Xy(</> 

.02169  71321 

.07429 

15920 

.27629 

23357 

1.1378 

38302 

-6 

.06315  81464 

.23428 

94293 

.96935 

98872 

4.5748 

10040 

-5 

.07820  10435 

.29577 

38500 

1.2557 

42634 

6.1234 

69370 

-4 

.09966  37988 

.38597 

79652 

1.6911 

30853 

8.5814 

32049 

-3 

.13059  80027 

.52048 

20797 

2.3690 

78533 

12.613 

56291 

-2 

.17548  64274 

.72384 

98638 

3.4491 

59992 

19.453 

18457 

-1 

.24077  73591 

1.0346 

12281 

5.2064 

25434 

31.446 

96885 

0 

.33548  43650 

1.5126 

74268 

8.1193 

82192 

53.177 

35907 

1 

.47177  85899 

2.2502 

13248 

13.028 

97303 

93.844 

08351 

2 

.66549  58331 

3.3878 

89404 

21.432 

66966 

172.47 

94233 

3 

.93653  73661 

5.1397 

44742 

36.043 

12414 

329.71 

05887 

4 

1.3092  61848 

7.8329 

42384 

61.872 

58912 

654.91 

96171 

5 

1.8130  90353 

11.972 

72241 

108.36 

34174 

1350.0 

16099 

6 

2.4835  95568 

18.348 

79645 

193.58 

95155 

2880.9 

67995 

7 

3.3642  83219 

28.207 

60768 

352.54 

26287 

6341.2 

44308 

8 

4.5091  48759 

43.527 

50601 

653.47 

59357 

14330. 

61359 

9 

5.9860  00189 

67.454 

60329 

1230.1 

76332 

33095. 

11689 

10 

7.8805  39970 

104.99 

11600 

2345.8 

02719 

77769. 

26471 
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APPENDIX  F.  INDEFINITE  INTEGRAL  OF  EXPONENTIAL  INTEGRAL 
The  integral  of  interest  here  is,  for  z  >  0, 

m  m 

®ll^*^  ■  J  ^  ®l^^^  “  J  ^  exp(-xt)  .  (F-1) 

X  1 

A  series  expansion  for  E^^(x) ,  which  is  good  for  ssiall  x,  is 
given  by  [9;  pages  258  -  259]: 


Eii(x)  -  i[ln(x)+Yl^  +  15  +  C  for  X  >  0  .  (F-2) 

However,  for  large  x,  this  series  loses  significance,  and  an 
asymptotic  expansion  is  required. 

We  begin  with  the  second  form  in  (F-1)  and  express 

•  2ZZ  (-1)*^'^  Yn  (t-1)"  for  |t-l|  <  1  ,  (F-3) 
n»l 

where 

Yjj*  **’  ■’’n  nil.  (F-4) 


Then,  as  X  +«, 


E 


W 

jj(x)  ~  J  dt  exp(-xt)  EZ  (-1)""^  Yn  (t-1)” 

■t  n"l 


exp(-x)  ^n 


x^  n^  (-x)""^ 


(F-5) 
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The  leading  terms  of  (F-5)  are  given  by 


Eii(x) 


as  X 


(F-6) 


A  program  for  for  x  >  0,  incorporating  series 

expansion  (F>2)  and  asymptotic  expansion  (F-5),  is  given  below  in 
BASIC.  The  switch  between  the  two  expansions  is  done  at  x  •  13, 
which  results  in  6  significant  digits  in  result  B^^tx).  A 
recursion  for  n!  Y||  is  built  into  the  program  in  lines  190  -  200. 


10  DBF  FNBll(X) 

20  IF  X>13.  THBN  150 
30  Tol«l.B>14 

40  T«1.0G(X)+. 57721566490153286 

50  S«.5*T*T+PI*PI/12. 

60  T-1. 

70  Sn«0. 

80  FOR  N«1  TO  100 
90  T--T*X/N 
100  R»T/(N*N) 

110  Sn<BSn+R 

120  IF  ABS(R)<Tol*ABS(Sn)  THEN  140 

130  NEXT  N 

140  RETURN  S+Sn 

150  T--1./X 

160  A-l.ElO 

170  H-G-P-S-1. 

180  FOR  N-2  TO  100 
190  H«H*(N-1) 

200  G«G*N-i-H 
210  P«P*T 
220  Ao«A 
230  A«G*P 

240  IF  ABS(A)>ABS(Ao)  THEN  270 

250  S-S+A 

260  NEXT  N 

270  S»S-.5*Ao 

280  RETURN  S*EXP(-X)*T*T 

290  FNEND 
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APPENDIX  G.  DIRECT  EVALUATION  OF  EXACT  EXCEEDANCE  DISTRIBUTION 


The  exact  characteristic  function  of  output  *  g(z^)  in 
(4),  of  the  nonlinear  device  in  (3),  can  be  found  with  the  aid  of 
(2).  For  the  case,  of  the  signal-present  bins  with  equal 
signal  powers  ^  *  S,  it  is  given  by  the  integral  expression 


exp(il 


exp(i5g(Xjj)) 


J  du  qj(u)  exp(i5g(u)) 


■  1  -  exp(-ax^)  +  J  du  a  exp(-au)  exp(iC[u  -  1  -  ln(u)])  ,  (G-1) 

*o 

where  a  ■  1/(1  *  S) •  For  case  Bq,  replace  a  by  1.  The 
characteristic  function  of  decision  variable  z  in  (4),  for  signal 
present  in  M  bins,  is  then  given  by 

This  result  can  be  evaluated  numerically  by  brute  force  from 
(G-1)  for  numerous  values  of  C,  and  used  in  the  fast  accurate 
procedure  of  [7]  to  directly  obtain  the  exceedance  distribution 
function  of  decision  variable  z,  at  least  for  some  cases  of 
breakpoint  x^  and  large  N.  An  example  for  S^O,  x^*!,  and 
N  >  1024  is  given  in  figure  G-1,  where  we  evaluated  integral 
(G-1)  at  increment  >  k/120,  up  to  C  ■  This  plot  is  the 
false  alarm  probability  and  should  be  compared  with  the  corres¬ 
ponding  simulation  cuirve  in  figure  3.  Agreement  is  excellent, 
but  we  are  now  eUdle  to  investigate  down  to  the  10  ^  level. 
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